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The  NPS  computer  simulation  model  was  used  to  investi- 
gate the  sputtering  by  1.0  keV  argon  bombardment  of  clean 
and  oxygen  reacted  monocrystalline  titanium,  vanadium  and 
niobium.   Variations  in  yield,  energy  and  angular  distri- 
butions and  in  particular,  multimer  formation,  were  studied 
as  a  function  of  oxygen  density  and  location.   Simulation 
results  show  a  significant  decline  in  substrate  yield  as 
the  oxygen  coverage  is  increased,  regardless  of  location. 
Further,  there  is  a  marked  preference  for  multimer  forma- 
tion by  lattice  fragmentation,  rather  than  by  recombination, 
for  all  three  metals.   The  percentages  of  multimers  formed 
by  fragmentation  were  found  to  increase  with  increasing  oxy- 
gen density  for  all  three  metals.   Previous  results  on  face- 
centered  cubic  crystals  concerning  the  yield  per  layer  and 
the  effects  of  channeling  were  confirmed  for  body-centered 
cubic  and  hexagonal  closed-packed  crystals.   For  compari- 
son, simulation  runs  were  repeated  at  ion  energies  of  0.5 
and  2.0  keV  for  a  titanium  with  oxygen  atoms  located  in 
three- fold  sites  with  a  C(2x2)  density. 
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I.   BACKGROUND 

A.   OVERVIEW 

Sputtering  is  the  ejection  of  atoms  from  a  target  by 
bombardment  with  an  energetic  projectile  ion.   The  incoming 
ion  collides  with  atoms  of  the  solid,  thereby  losing  a 
portion  of  its  energy  and  momentum.   If  the  energy  gained 
by  a  struck  atom  exceeds  the  binding  energy  holding  it  in 
position  in  the  lattice  a  primary  recoil  atom  will  be 
created.   This  recoil  atom  will  then  collide  with  other  atoms 
in  the  material,  distributing  its  energy  via  a  collision 
cascade.   An  atom  at  or  near  the  surface  will  be  sputtered 
if  the  energy  associated  with  the  normal  component  of  its 
velocity  exceeds  the  surface  binding  energy  of  the  material. 

The  phenomenon  of  sputtering  can  be  subdivided  into  two 
principle  categories.   Transmission  sputtering  is  the  ejec- 
tion of  atoms  from  the  rear  of  a  thin  target  following  pene- 
tration by  the  ion.   Sufficient  energy  must  be  transported 
through  the  target  to  allow  atoms  to  overcome  their  binding 
energy  and  escape.   Back-sputtering  is  the  ejection  of  atoms 
from  the  front  of  the  target  and  depends  upon  sufficient 
energy  being  deposited  in  the  surface  layers  to  allow 
ejection.   This  thesis  will  be  concerned  with  the  latter. 

First  discovered  by  Grove  in  1853  and  Faraday  in  1854 
as  the  deposition  of  metal  atoms  on  the  glass  walls  of  a 
gas  discharge  tube,  it  was  not  until  19  02  that  Goldstein 


[Ref.  1]  presented  evidence  that  the  sputtering  effect  was 
caused  by  positive  ions  of  the  discharge  impacting  on  the 
metal  cathode.   Since  its  principle  effect  was  to  contaminate 
experimental  environments  and  to  erode  apparatus,  the  majority 
of  early  work  centered  around  ways  to  control  or  eliminate 
what  was  considered  to  be  an  undesirable  side  effect.   The 
necessity  of  increasing  the  life  expectancies  of  high  voltage 
vacuum  tubes  and  the  recognition  of  the  possible  applications 
of  sputtering  to  thin  film  coatings  gradually  spurred  inter- 
est in  understanding  the  mechanics  involved. 

Early  theoretical  work  was  hampered  by  a  paucity  of 
experimental  data.   Reproducibility  of  results  was  also  a 
major  problem  due  to  an  incomplete  understanding  of  the 
factors  involved.   Since  the  sputtering  yield,  Y  (defined  as 
the  number  of  atoms  ejected  per  incident  ion) ,  is  critically 
sensitive  to  ambient  pressure  and  surface  contamination,  as 
well  as  the  flux  density  and  angle  of  incidence  of  the  beam, 
progress  was  dependent  on  the  development  and  refinement  of 
experimental  techniques . 

B.   EARLY  EXPERIMENTAL  WORK 

The  effect  of  pressure  on  yield  was  first  demonstrated 
by  Penning  and  Moubis  [Ref.  2] .   They  found  that  an  increase 
in  pressure  increased  the  collision  frequency  between  the 
ejected  atoms  and  surrounding  gases,  resulting  in  the  back- 
scattering  of  the  escaping  atoms  back  onto  the  target 

-5 
surface.   By  keeping  background  pressure  below  10    Torr, 
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they  were  able  to  obtain  reproducible  results  for  ion  ener- 
gies in  excess  of  500  eV.   Below  500  eV,  the  scouring  action 
of  the  incident  beam  was  insufficient  to  maintain  a  clean 
surface. 

The  elimination  of  adsorbed  gases  and  surface  oxide 
layers  as  a  critical  factor  in  determining  the  true  yield 
was  described  by  Arifov,  et  al.,  in  1963  [Ref.  3].   Yonts 
and  Harrison  also  presented  evidence  [Ref.  4]  that  surface 
recontamination  from  background  gases  was  a  significant 
factor  in  quantitative  sputtering  yield  measurements. 
Criteria  were  therefore  developed  for  conducting  'clean1 
experiments . 

In  addition  to  pressure  and  surface  cleanliness,  the 
following  conditions  must  be  met  in  order  to  obtain  repro- 
ducible yields: 

(1)   The  ions  must  strike  the  target  at  a  known  angle 
of  incidence. 

C2)   The  energy  spread  of  the  beam  must  be  small. 

(3)  All  ions  must  be  uniformly  charged  and  mass  separated, 

(4)  On  monocrystalline  surfaces  the  lattice  orientation 
must  be  adequately  described. 

C.   THEORETICAL  DEVELOPMENT 

Pioneering  work  on  the  concept  of  an  individual  sputter- 
ing event  on  an  atomic  scale  was  conducted  by  Stark  [Ref.  5] 
Applying  the  conservation  laws  governing  elastic  collisions, 
he  proposed  two  theories  to  explain  known  experimental 


11 


results.   One,  the  "Hot-Spot"  model,  considered  the  sputtered 
atoms  to  be  the  result  of  evaporation  of  target  material 
from  an  atomically  small  region  due  to  ion  induced  local 
heating.   In  his  "Collision"  model,  sputtering  events  were 
the  result  of  a  series  of  binary  collisions  initiated  by 
a  single  ion.   His  models  were  modestly  successful  in  pre- 
dicting the  energy  distribution  of  the  sputtered  particles. 
Further  refinements  on  the  hot-spot  model  were  attempted  by 
von  Hippie  [Ref .  6] ,  but  only  after  the  demonstration  by 
Wehner  [Ref.  7]  of  the  effects  of  crystal  structure  on  the 
yield  did  it  become  apparent  that  local  heating  alone  could 
not  explain  the  mechanisms  involved  in  sputtering. 

Wehner' s  observation  of  spot  patterns  (angular  distri- 
bution patterns  characteristic  of  the  crystallographic  proper- 
ties of  the  substrate)  obtained  from  sputtering  of  mono- 
crystalline  targets  revived  interest  in  the  collision  approach 
In  19  52  Keywell  [Ref.  8]  made  a  first  attempt  to  formulate 
Stark's  collision  model  in  terms  of  the  neutron  transport 
model  developed  for  nuclear  work.   Harrison  then  developed 
a  theory  based  on  probability  concepts  as  expressed  by  the 
idea  of  collision  cross-sections  [Ref.  9 J .   However,  due  to 
the  large  number  of  unknown  parameters,  these  theories  were 
limited  in  their  quantitative  application. 

As  a  further  attempt  to  account  for  the  angular  distri- 
bution of  ejecta,  Silsbee  [Ref.  10 J  proposed  a  focused 
collision  model  which  allowed  the  transport  of  momentum  in 
crystals  along  preferred  directions.   Available  experimental 
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results  showed  that  the  yield  from  monocrystals  depended 
sensitively  on  the  crystallographic  orientation  of  the  inci- 
dent beam.   This  was  explained  in  terms  of  the  "holes"  seen 
in  crystal  models  when  viewed  from  different  planes  and 
led  to  the  development  of  various  "Transparency"  models. 
While  focusing  has  been  shown  to  be  a  contributing  factor 
in  cascade  development  within  a  crystal,  it  is  inadequate 
to  explain  Wehner  spots. 

As  an  alternative  to  the  focusing  collision  theory, 
Lehmann  and  Sigmund  [Ref .  10]  proposed  a  model  based  on  the 
Boltzman  transport  equations  that  required  the  target  surface 
to  have  an  ordered  structure,  but  not  long  straight  rows  of 
atoms  intersecting  the  surface  as  in  the  focusing  model. 
Emphasis  is  placed  on  the  role  of  the  surface  structure  and 
in  particular,  surface  binding  energies  in  monocrystalline 
sputtering.   As  an  ion  strikes  the  surface  a  collision  cas- 
cade is  created  via  binary  interactions.   An  atom  in  the 
cascade  will  sputter  if  its  momentum  component  normal  to 
the  surface  has  sufficient  energy  associated  with  it  to 
overcome  the  surface  potential  barrier.   Thompson  {Ref.  12] 
expanded  on  this  idea  and  proposed  that  the  surface  attrac- 
tion for  the  escaping  particle  causes  a  refraction  of  its 
velocity  vector  away  from  the  normal,  resulting  in  a  dis- 
torted angular  distribution  of  ejected  particles.   Signifi- 
cant agreement  with  experimental  results  was  found  when 
projectile  energies  were  much  greater  than  the  surface  bind- 
ing energy.   This  approach  also  provided  correct  information 
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concerning  both  the  angular  and  the  energy  distributions 
from  polycrystalline  targets.   Thompson's  [Ref.  12]  work 
also  incorporated  the  idea  of  focusons  (momentum  transport 
without  mass  displacement)  into  the  theory.   Robinson  [Ref. 
13]  has  taken  a  similar  approach  while  accounting  for 
channeling. 

While  these  theories  have  provided  mathematical  tools 
to  aid  experimentalists,  they  have  been  rather  unsuccessful 
in  predicting  yields  at  low  projectile  energies  and  from 
single  crystal  targets.   Binary  collision  models  have  thus 
far  been  unable  to  account  for  the  predominance  of  sputtered 
atoms  from  surface  layers  of  the  target. 

In  general,  these  theories  have  attempted  to  reduce  what 
is  essentially  a  many-body,  multiple  interaction  process 
into  forms  that  are  analytically  tractable.   An  assortment 
of  refinements  followed,  attempting  to  explain  observed 
behavior.   Full  scale  attack  on  the  multiple  interaction 
problem  had  to  await  the  arrival  of  high  speed  computers  with 
sufficient  storage  to  handle  the  hundreds  or  thousands  of 
atoms  in  a  cascade. 

D.   COMPUTER  SIMULATION  MODELS 

Paralleling  the  theoretical  and  experimental  progress  in 
the  understanding  of  sputtering  was  the  effort  to  model  the 
basic  mechanisms  through  the  use  of  computer  simulation.   In 
1960  Gibson,  Goland,  Milgram  and  Vineyard  [Ref.  14J  developed 
a  computer  model  to  simulate  the  motion  of  atoms  in  a  copper 
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monocrystal  following  the  impact  of  an  ion  on  a  single 
target  atom.   Atomic  interactions  were  treated  as  binary 
collisions  and  the  resulting  cascades  were  developed  by 
applying  Newton's  equations  of  motion.   Ion  penetration 
studies  by  Robinson  and  Oen  [Ref .  15]  using  a  similar 
program  led  to  the  discovery  of  channeling  which  was  later 
experimentally  verified. 

In  19  67,  based  on  earlier  preliminary  studies,  Harrison, 
Levy,  Johnson  and  Effron  [Ref.  16]  used  a  computer  to  simu- 
late the  bombardment  of  a  single  copper  crystal  by  argon. 
An  important  finding  of  this  work  was  that  the  yield  was 
almost  entirely  confined  to  the  first  three  layers  of  the 
top  surface.   Subsequent  refinements  of  the  model  [Refs.  17- 
20]  included  the  addition  of  potential  functions  that  allowed 
for  attractive,  as  well  as  repulsive  interactions,  creating 
a  dynamically  stable  crystal.   Harrison's  results  showed  that 
focusons  were  not  a  major  contributor  to  the  yield. 

In  19  78  Garrison,  Winograd  and  Harrison  [Ref.  21]  pub- 
lished the  results  of  a  comprehensive  study  of  atomic  and 
molecular  ejection  from  a  copper  crystal  with  adsorbed  oxy- 
gen atoms .   Their  simulation  of  the  system  provided  a  de- 
tailed picture  of  the  mechanisms  of  ejection  and  enabled 
them  to  determine  whether  molecules  were  ejected  as  a  cluster 
from  the  surface  or  combined  together  in  the  space  above  the 
target.   They  were  also  able  to  determine  the  effects  of 
adatom  placement  on  molecule  formation.   Such  results,  when 
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coupled  with  laboratory  data,  should  lead  to  a  better  under- 
standing of  actual  adatom  locations  on  chemically  reacted 
surfaces. 

While  much  work  remains  to  be  done,  particularly  with 
respect  to  constructing  interatomic  potentials  that  accurately 
represent  the  physical  processes  involved,  computer  simula- 
tion permits  the  study  of  sputtering  phenomena  with  a  preci- 
sion and  flexibility  that  is  impossible  to  duplicate  in  the 
laboratory.   While  it  cannot  replace  actual  experimentation, 
it  can  lead  to  a  better  understanding  of  experimental  re- 
sults and  provide  a  testing  ground  for  theoretical  concepts . 
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II.   OBJECTIVES 

The  primary  objective  of  this  thesis  is  to  study  the 
sputtering  of  titanium,  vanadium  and  niobium  monocrystals 
using  the  NPS  computer  simulation  model,  QRAD.   In  addition, 
changes  in  the  yield  and  multimer  formation  are  studied  when 
the  surface  is  reacted  with  oxygen. 

Yu  [Refs.  22,23]  has  conducted  comprehensive  research 
on  clean  and  oxygen  reacted  titanium,  niobium  and  vanadium 
in  an  attempt  to  link  several  earlier  theories  on  the 
mechanisms  responsible  for  ionic  multimer  formation.   In 
the  Lattice  Fragmentation  Model  first  proposed  by  Benning- 
hoved  [Ref.  24]  molecular  species  Cmultimers)  are  emitted  as 
complete  fragments  of  the  surface  lattice  during  ion  bombard- 
ment.  In  the  Recombination  Model  [Refs.  18,25]  only  single 
atoms  are  ejected.   Molecules  are  then  formed  through  the 
attractive  interaction  between  sputtered  atoms  from  the 
same  sputtering  event,  establishing  their  identities  above 
the  surface.   The  charge  state  of  the  molecules  is  not  taken 
into  account  in  these  models. 

Yu  has  reported  experimental  evidence  [Ref.  23]  that  the 
emission  mechanism  for  a  sputtered  molecular  species  may  be 
a  function  of  its  charge  state.   He  proposed  that  for  low 
oxygen  coverages  (<  one  monolayer)  the  formation  of  positively 
charged  titanium  and  niobium  molecules  (MO  ,  M0-)  favored 
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the  recombination  model,  while  the  formation  of  negatively 
charged  species  (M0~,  M0~  and  M0~)  followed  the  fragmenta- 
tion process  in  all  three  metals.   In  the  case  of  vanadium 
it  appeared  that  the  positive  oxide  species  were  not  formed 
exclusively  by  recombination,  but  by  a  combination  of  the 
two  methods. 

Formation  of  positively  charged  ions  was  also  found  to 
be  strongly  dependent  on  oxygen  coverage.   When  the  number 
of  oxygen  atoms  neighboring  a  metal  atom  increases,  the 
frequency  of  sputtering  an  M+  and  0   adjacent  to  each  other 
naturally  increases.   Consequently  the  percentage  of  M0+ 
produced  by  fragmentation  will  increase. 

The  formation  of  negative  molecular  ions  by  lattice 
fragmentation  was  found  to  be  insensitive  to  oxygen  coverage 
due  to  the  kinetics  involved  in  the  formation  proces .   A 
molecule  formed  by  recombination  is  usually  left  with  excess 
energy.   If  this  energy  approaches  the  electron  affinity 
of  the  molecule  the  negative  ion  will  attempt  to  stabilize 
itself  by  ejecting  an  electron.   If  the  excess  energy  is 
much  greater  than  this,  one  approaches  the  dissociation 
energy,  and  the  molecule  will  be  prevented  from  forming  at 
all.   Thus  the  formation  of  a  stable  negative  ion  is  most 
favorable  when  the  components  are  close  neighbors  on  the 
surface  and  eject  with  small  relative  velocities  as  in  the 
lattice  fragmentation  process.   Increasing  the  oxygen  cover- 
age will  only  increase  the  probability  of  forming  higher 
order  oxide  multimers . 
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Some  discrepancies  were  noted  in  Yu's  results  for  posi- 
tive ions  and  were  attributed  to  the  recombination  of  metal 
atoms  and  oxygen  from  adjacent  sites.   Thus  the  results  were 
indistinguishable  from  fragmentation  at  higher  oxygen 
coverages.   All  nine  of  the  negative  molecular  ion  species 
investigated  by  Yu  showed  strong  evidence  of  formation  by 
fragmentation.   Since  the  energy  balance  for  recombined 
negatively  charged  ions  indicates  that  they  tend  to  neutralize 
by  electron  ejection,  this  leaves  fragmentation  as  the 
principal  source  of  molecule  formation. 

It  should  be  noted  that  Yu's  experimental  procedure 
utilized  an  ion  beam  with  a  45  degree  angle  of  incidence  to 
the  polycrystalline  surface.   The  mass  spectrometer  detector 
was  also  positioned  at  45  degrees  with  respect  to  the  surface, 
at  right  angles  to  the  ion  beam.   Thus  only  a  portion  of  the 
total  yield  was  sampled.   The  data  from  the  simulation  studies 
that  follow  are  based  on  the  total  yield  from  ion  bombardment 
perpendicular  to  monocrystal  surface.   These  data  include 
neutral  atoms,  as  well  as  ions,  which  are  not  collected  by 
the  mass  spectrometer. 

While  determination  of  the  charge  state  of  molecules  is 
not  possible  with  the  computer  model,  the  various  mechanisms 
of  ejection  can  be  studied  in  detail  to  determine  whether 
recombination  or  fragmentation  is  the  major  contributor  to 
the  multimer  yield.   Ejection  times  of  the  various  constitu- 
ents are  also  available  for  correlation. 
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Additionally  the  energy  and  angular  distributions  of 
ejected  particles,  which  can  be  related  to  experimental  SIMS 
data,  will  be  studied.   Since  the  actual  sites  of  adatom 
adsorption  have  yet  to  be  determined,  the  effects  of  oxygen 
coverage  and  location  on  yield  and  multimer  formation  will 
also  be  investigated. 
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III.   MODEL  DEVELOPMENT  AND  SIMULATION 

A.   COMPUTER  MODEL 

QRAD  is  a  computer  simulation  program  designed  to  track 
the  motions  of  atoms  in  a  microcrystal  following  ion  bombard- 
ment.  Trajectories  are  calculated  by  solving  Hamilton's 
equations  of  motion  using  an  average  force  method  [Ref.  26]. 
It  is  a  multiple  interaction  (MI)  simulation  following  time 
step  logic.   Once  the  atomic  masses,  crystal  characteristics, 
interatomic  potentials  and  ion  energy  have  been  initialized, 
the  program  will  develop  and  follow  the  resulting  collision 
cascade  until  the  energy  of  the  atoms  reaches  a  predetermined 
minimum.   (See  reference  24  for  a  comprehensive  discussion 
of  sputtering  simulations . ) 

Initially  an  ion  of  given  mass  is  projected  toward  a 
surface  impact  point  with  a  specified  energy  and  angle  of 
incidence.   As  the  cascade  develops,  the  positions  and  veloci- 
ties of  the  ion  and  lattice  atoms  are  calculated  at  the  end 
of  each  timestep.   A  timestep  is  determined  by  the  velocity 
of  the  most  energetic  particle  so  that  its  displacement  does 
not  exceed  0  J.  lattice  unit  CLU) .    To  expedite  calculations, 
particles  are  classified  as  either  moving  or  non-moving  de- 
pending on  their  energies  and  the  forces  exerted  on  them 


The  lattice  unit  is  a  convenient  method  of  expressing 
distances  when  dealing  with  different  crystal  types.   An 
LU  is  equal  to  the  lattice  constant  'a*  divided  by  root  two 
for  hep  crystals  and  is  equal  to  a/2  for  bec  crystals. 
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by  neighbor  atoms.   Atoms  which  rise  above  the  surface  with 
sufficient  momentum  are  reclassified  and  assigned  to  a 
tentative  list  of  ejected  atoms.   A  trajectory  is  terminated 
when  the  energy  of  the  most  energetic  atom  is  such  that  the 
probability  that  either  it,  or  other  atoms,  will  escape  the 
surface  is  very  small.   Final  positions  and  velocities  are 
then  recorded.   Multiple  trajectories,  simulating  bombard- 
ment by  an  ion  beam,  can  be  run  by  computing  trajectories 
for  a  set  of  impact  points  chosen  to  sample  a  representative 
area  of  the  crystal  surface  as  a  function  of  its  symmetry. 
Data  on  the  sputtered  atoms  is  then  analyzed  by  separate 
programs  (AN83  and  PLOTAN)  to  determine  multimer  formation, 
ejection  angle  information,  etc. 

As  mentioned  previously,  multimer s  may  result  either 
from  lattice  fragmentation  or  by  recombination  above  the 
surface.   A  multimer  is  considered  formed  when  its  constitu- 
ent atoms  are  in  a  favorable  spatial  position  with  suffi- 
ciently low  relative  kinetic  energy  to  permit  bonding.   The 
total  energy  of  the  resulting  system  must  be  negative  and 
is  an  indication  of  the  group's  binding  stability.   Absolute 
multimer  yield  is  strongly  dependent  on  the  pair-potential 
in  use  and  the  range  specified  for  atomic  interactions 
(Rp) ,  however  the  mechanisms  for  cluster  formation  have  been 
found  to  be  insensitive  to  the  choice  of  the  potential 
function  [Ref .  27]  . 

Global  results,  suitable  for  comparison  with  experimental 
data,  are  obtained  by  averaging  the  yields  from  each  run 
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over  a  number  of  trajectories.   While  the  number  of  trajec- 
tories is  quite  small  (100-200,  depending  on  symmetry 
considerations) ,  the  results  have  been  found  to  vary  by 
less  than  10%  when  compared  with  other  independently  chosen 
sets  of  impact  points  [Ref .  2  8] .   These  global  results  are 
used  in  determining  the  number  of  atoms  ejected  per  single 
ion  (ASI  yield) ,  energy  and  angular  distributions,  as  well 
as  the  spatial  distribution  of  yield  from  the  impact  zone. 

The  remainder  of  this  chapter  is  devoted  to  the  develop- 
ment of  the  potential  functions  and  other  parameters  used 
by  QRAD.   Adatom  nomenclature  and  placement  are  discussed 
in  Section  D.   Trajectory  sets  and  crystal  impact  points  are 
covered  in  Sections  F  and  G.   The  chemical  properties  of 
titanium,  vanadium  and  niobium  are  covered  in  the  following 
sections. 

B.   TITANIUM 

Titanium  is  a  dimorphic  transition  metal  ranking  ninth 
in  abundance  of  the  elements  making  up  the  earth's  crust. 
Below  about  880 °C  it  exists  as  a  silvery  gray  hexagonal 
close-packed  (hep)  crystalline  material,  while  above  880 °C 
it  undergoes  a  phase  transformation  to  a  face-centered 
cubic  (fee)  structure.   The  important  physical  and  chemical 
properties  of  titanium  are  listed  in  Table  1. 

Titanium  reacts  with  oxygen  to  produce  four  well  defined 
oxides,  the  monoxide  TiO,  the  sesquioxide  Ti203,  the  dioxide 
or  titanic  acid  Ti02  and  the  trioxide  or  pertitanic  acid 
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Ti03*   Ti2°5'  Ti3°5  and  Ti7°l2  have  also  been  reported,  but 
may  in  fact  be  metastable  combinations  of  the  others  [Ref. 
29]  . 

At  least  two  distinct  forms  of  oxygen  coverage  have  been 
observed  on  Ti(OOl)  surfaces. 

1 .  Type  1  begins  to  form  at  low  oxygen  exposures  of 

2 

approximately  1.0  L  and  is  characterized  by  a  low  work 

function,  very  low  electron  stimulated  desorption  (ESD) 
and  photon  stimulated  desorption  (PSD)  yields  and  a  surface 
binding  energy  of  5.9  eV  as  determined  by  UV  photoemission 
spectroscopy  (UPS).   Plateau,  et  al.,  had  previously  deter- 
mined this  value  to  be  5.7  ±  0.2  eV,  also  using  UPS  and 
ESCA  [Ref.  38] . 

2.  As  the  oxygen  exposure  is  increased  (-50  L)  a  second 
type  of  oxygen  coverage  is  observed  which  has  a  higher  work 
function,  high  PSD  yield  and  additional  binding  energy  in 
UPS.   The  effects  of  oxygen  adsorption  appear  to  saturate 
at  about  70  L. 

Hanson  [Ref.  37]  concluded  that  type  1  may  correspond 
to  oxygen  adatoms  lying  below  the  surface  plane  or  in  high 
coordination  sites,  possibly  the  three-fold  hollows,  which 
would  account  for  the  low  desorption  rate. 

Jonker  [Ref.  39]  has  also  observed  two  states  of  oxygen 
on  the  surface.   The  first  is  detected  at  low  coverages  and 
produces  a  well-ordered  P(2x2)  low  energy  diffraction 


0  -6  -4 

1.0  L  (Langmuir)  =  10    Torr-sec  =  1.33  x 10    Pa-sec. 
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(LEED)  pattern.   At  higher  oxygen  exposures  a  second  state 
emerges  which  has  a  disordered  appearance.   This  was  ex- 
plained by  assuming  that  the  sticking  probability  coefficient 
for  oxygen  is  relatively  constant  up  to  a  quarter  monolayer 
coverage.   At  this  point  interference  from  previously 
adsorbed  oxygens  would  cause  it  to  change  abruptly.   The 
sticking  coefficient  would  then  remain  constant  for  increas- 
ing exposures  until  the  density  of  oxygen  favored  a  change 
in  the  chemical  bonding,  resulting  in  the  formation  of  a 
more  stable  metal  oxide  in  the  upper  layers  of  the  crystal. 

Kawasaki,  et  al.  [Ref.  40],  have  also  reached  the  conclu- 
sion that  oxygen  is  first  adsorbed  on  the  surface  at  low 
doses,  but  as  the  exposure  is  increased,  it  is  slowly  ad- 
sorbed into  the  surface.   This  corresponds  to  the  formation 
of  an  oxide,  and  is  accompanied  by  changes  in  the  crystal 
structure.   Their  results  also  indicate  that  at  the  lower 
exposures  oxygen  occupies  either  a  doubly  bonded  A-top  or 
a  two-fold  bridge  position.   The  A-top  position  is  also  sup- 
ported by  PSD  experiments  [Ref.  41]  in  which  the  oxygen 
atoms  were  observed  to  desorb  normal  to  the  surface,  imply- 
ing that  the  Ti-0  bonds  are  perpendicular  to  the  surface. 

Experimental  work  by  Singh,  et  al .  [Ref.  43],  confirms 
the  formation  of  an  oxide  as  saturation  is  reached,  with 
oxygen  atoms  incorporated  extensively  into  the  metal  lattice. 
(Here  saturation  is  defined  as  the  point  where  increasing 
the  oxygen  exposure  no  longer  has  an  effect  on  the  metal 
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crystal,  i.e.,  the  sticking  probability  is  equal  to  zero.) 
Since  the  sputtering  yield  comes  mainly  from  the  top  three 
layers,  modeling  of  an  oxygen  saturated  surface  would  have 
to  be  accomplished  by  using  a  titanium  oxide  crystal  rather 
than  titanium  with  oxygen  adatoms. 

In  view  of  the  uncertainty  of  the  actual  locations  of  the 
oxygen  atoms,  a  complete  set  of  trajectories  was  run  for  the 
A- top,  two-fold  (bridge)  and  three-fold  locations,  in  addi- 
tion to  those  on  the  clean  surface.   Figures  1  through  8 
present  the  various  adatom  locations  on  the  metal  crystal 
surface.   (See  Section  D  for  a  description  of  adatom  notation.) 

C.   VANADIUM  AND  NIOBIUM 

Vanadium  and  niobium  are  also  transition  metals,  but 
unlike  titanium,  they  have  a  body-centered  cubic  (bcc) 
structure.   Their  (110)  planes,  while  not  as  close  packed, 
do  however  have  a  hexagonal  arrangement  of  atoms  similar  to 
the  titanium  (001)  plane  (Fig.  9) .   Their  physical  and 
chemical  properties  (Tables  2  and  3)  are  also  quite  similar 
to  titanium,  hence  their  sputtering  mechanisms  should  bear 
some  correlation  to  that  of  titanium. 

Haas,  Jackson  and  Hooker  [Ref.  44 J  in  their  LEED  experi- 
ments on  niobium,  tantalum  and  vanadium  have  found  that  oxygen 
adsorption  is  the  same  on  all  three  metals  and  is  a  function 
of  oxygen  exposure.   At  low  oxygen  coverages  the  oxygen 
patterns  clearly  indicate  a  six-fold  symmetry  characteris- 
tic of  the  substrate  lattice.   But,  as  in  the  case  of 
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Fig.  1.  Titanium  (001)  surface  with  adatoms  in  a  P(2x2)  A-top 

configuration.  Rectangles  in  this  and  subsequent  figures 
indicate  impact  zones  used  in  the  simulation  runs. 
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Fig.   2.     Titanium  with  adatoms   in  a  C(2x2)  A-top  configuration 
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Fig.  3.  Titanium  with  adatoms  in  P(2x2)  bridge  configuration 
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Fig.  4.  Titanium  with  adatoms  in  C(2x2)  bridge  configuration 
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Fig.   5.     Titanium  with  adatoms  in  P(2x2)   three-fold  "A"  position 
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Fig.  6.  Titanium  with  adatoms  in  C(2x2)  three- fold  "A"  position 
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Fig.   7.     Titanium  with  adatoms  in  P(2x2)  three-fold  "B"  position 
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Fig.  8.  Titanium  with  adatoms  in  C(2x2)  three-fold  "B"  position 
dots  indicate  centers  of  atoms  in  the  second  layer. 


Small 
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Fig.  9.  Vanadium  and  Niobium  clean  (110)  surface 
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Fig.  10.  Vanadium  and  niobium  with  adatoms  in  P(2x2)  A-top  position 
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Pig.   11.     Vanadium  and  niobium  with  adatoms  in  C(2x2)  A-top  position 
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TABLE  2 

Physical  Properties  of  Vanadium  and  Parameters 
of  Associated  Potential  Functions 


Atomic  Number 

Atomic  Weight 

Melting  Point 

Boiling  Point 

Valance 

Crystal  Type 
a 

Lattice  Unit 

Re 

Cohesive  Energy 

0 

Gas  Phase  Bond  Lengths  (A) 

V-V 
Ti-0 
0-0 
Ar-0 

Gas  Phase  Bond  Strengths  (eV) 

V-V 
V-0 
0-0 
Ar-0 

Surface  Binding  Energy 

V-0 
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REF 

50.9415 

30 

1890  ±  10°  C 

30 

3380°  C 

30 

+2,3,4,5 

30 

bcc 

31 

3.0399  ±  0.0003  A 

1.51995  A 

Calc. 

2.6326  A 

31 

5.31  eV 

33 

2.6224 

30 

1.56  ±  .04 

34 

1.20752 

35 

3.31 

35 

2.511 

30 

6.462 

30 

5.1545 

35 

0.0049 

30 

5.52  eV 

Calculated, 

see  text. 
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TABLE  3 

Physical  Properties  of  Niobium  and  Parameters 
of  Associated  Potential  Functions 


Atomic  Number  41  REF 

Atomic  Weight  92.9064  30 

Melting  Point  2468  ±10°  C  30 

Boiling  Point  4742°  C  30 

Valance  +2,3,4(7) ,5  30 

Crystal  Type  bcc              0   31 

a  3.30656  ±  0.00002  A 

Lattic  Unit  1.65328  A  Calc . 

Re  2.8636  A  Calc. 

Cohesive  Energy  7.57  eV  33 

o 

Gas  Phase  Bond  Lengths  (A) 

Nb-Nb  2.8584  30 

Nb-0  1.69  34 

0-0  1.20752  35 

Ar-0  3.31  35 

Gas  Phase  Bond  Strengths  (eV) 

Nb-Nb  2.83  34 

Nb-0  7.805  35 

0-0  5.1545  35 

Ar-0  0.0049  30 

Surface  Binding  Energy 

Nb-0  6.67  eV  Calculated, 

see  text. 
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titanium,  the  reference  index  for  the  adatom  mesh  is  uncer- 
tain.  Two  locations  are  possible  on  the  vanadium  and  niobium 

(110)  faces;  the  A-top  and  a  rectangular  four-fold  site  not 
present  in  titanium.   In  the  four-fold  position  the  adatoms 
could  be  placed  either  equidistant  from  the  two  pairs  of 
neighbor  atoms,  or  off  to  one  side  in  a  three-fold  'plus' 
position.   Both  would  result  in  a  high  degree  of  bonding 
asymmetry.   Since  there  is  no  corresponding  site  on  the  Ti 

(001)  surface,  the  A-top  position  was  chosen  for  comparison 
with  titanium. 

No  information  is  available  on  the  binding  energies  of 
chemisorbed  oxygen  on  vanadium  or  niobium.   Since  all  three 
metals  possess  similar  chemical  properties,  surface  binding 
energies  were  calculated  by  scaling  the  gas  phase  M-0 
binding  energies  to  the  known  SBE  of  oxygen  on  titanium. 
This  resulted  in  estimated  SBE's  of  5.52  and  6.67  eV 
respectively . 

D.   ADATOM  COVERAGE 

Computations  were  run  at  two  specific  coverages,  P(2x2) 
and  C(2x2),  for  each  of  the  adatom  locations.   These  cover- 
ages correspond  roughly  to  a  quarter  and  a  half  monolayer 
of  oxygen  respectively. 

The  designations  for  surface  coverages  were  standardized 
in  an  article  by  E.A.  Wood  in  1964  [Ref.  45].   In  this  nota- 
tion the  three  dimensional  unit  cell  is  reduced  to  a  diperiodic 
structure  referred  to  as  a  net  in  which  the  unit  area  is  the 
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unit  mesh.   The  fourteen  Bravais  space  lattices  are  thus 
reduced  to  five  basic  nets  pictured  below. 


Hexagonal 
y   =  120',  a  =  b 


Oblique 

Y  /90',  a  /  b 


Rectangle 
Y  =  90',  a  f   b 


Centered  Rectangle 
Y  =  90*,  a  f   b 


The  letter  "P"  is  used  to  denote  the  primative  mesh  with 
equivalent  points  at  the  four  corners,  but  none  within  the 
mesh  itself.   The  letter  "C"  is  used  to  identify  a  unit  mesh 
that  contains  a  central  point.   Numerals  are  used  to  indi- 
cate the  repetition  interval  in  the  X  and  Z  directions, 
referenced  to  the  substrate  unit  mesh.   Meshes  rotated  from 
the  substrate  principal  directions  are  designated  as  above 
followed  by  the  angle  of  rotation.   NaCl (100)  PC5xl)-45 
would  represent  a  primative  mesh  located  on  the  (100). 
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surface  of  NaCl,  five  units  by  one  unit  in  area  and  oriented 
in  the  <110>  direction.   Meshes  formed  on  other  than  sub- 
strate atoms  are  further  identified  by  their  origin  point, 
two-fold  bridge,  etc.   (For  simplicity  the  four  number 
designation  normally  used  for  hep  crystals  has  been  abbre- 
viated to  three  digits;  i.e.,  (0001)  =  (001).   The  numbers 
themselves  refer  to  the  principal  crystal  directions  (hkil) 
where  i  =  h  +  k.) 

In  the  case  of  the  titanium  (001)  surface  and  the  (110) 
surfaces  of  vanadium  and  niobium,  the  unit  mesh  is  a  120 
degree  rhombus.   The  unit  mesh  for  the  A-top  positions  is 
centered  above  the  substrate  atoms  and  is  displaced  for  the 
two  and  three-fold  positions.   Figures  1-11  depict  the  various 
coverages  used,  with  the  adatoms  scaled  according  to  their 
covalent  radii.   The  two-fold  locations  are  abbreviated  "BR" 
(bridge)  and  the  three-fold  sites  by  "TF".   Since  there  are 
two  different  three-fold  sites  on  the  hep  (001)  surface,  the 
additional  identifier  "TFA"  or  "TFB"  is  used  to  denote 
whether  the  adatom  is  fixed  "above"  an  atom  in  the  second 
layer  of  the  substrate  or  centered  over  a  channel. 

E.   POTENTIAL  FUNCTIONS 
1.   General 

A  variety  of  potential  function  forms  are  available 
in  this  program  to  describe  the  interatomic  forces.   The 
merits  of  each  have  been  hotly  debated  in  the  literature, 
but  Harrison  has  found  [Ref.  42]  that  the  global  results  are 
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not  strongly  dependent  on  the  particular  form  in  use.   Indi- 
vidual parameters  may  be  adjusted,  if  necessary,  so  that 
absolute  values  match  experimental  results,  but  the  basic 
mechanisms  involved  remain  unchanged.   The  four  potential 
functions  used  are  as  follows: 

a.  Born-Mayer ,  which  is  strictly  repulsive  having  the 
form 

V(r)   =  A  exp(-br)  ; 

b.  Moliere,  which  is  an  approximation  to  the  Thomas-Fermi 
screening  function,  and  is  also  repulsive 


V(r)   =   (Z  Ze2/r) {0.35  exp  (-0.3r/a)  +  0  .  55  exp  (-1 .2r/a) 

1  « 


+  0.1  exp  (-6.  Or/a)  }  , 


where: 


=   0.8853aQ/(z:[/2  +  Z^/2)  2/3    (Firsov  radius) 


and 


a   =   h2/4ir2me2   =   0.5292  A.     (Bohr  radius) 


The  parameter  "a"  may  be  varied  to  produce  a  "modified" 
Moliere  function. 
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c.  Morse,  which  is  both  repulsive  and  attractive  having 
the  form 

V(r)   =   D  exp  {-2a(r-r  )  }  -  2D  exp{-a(r-r  )}  , 

where  rQ   is  the  equilibrium  separation  distance  between 
atoms,  De  is  the  well  depth  and  alpha  is  used  to  scale  the 
distance  in  the  potential  function;  and 

d.  Composite  Moliere-Morse,  which  consists  of  the  two 
potential  functions  joined  together  smoothly  by  a  cubic 
spline. 

The  adjustable  parameters  can  be  related  to  the  ion 
"size"  and  a  collision  "hardness",  which  in  the  case  of  the 
Born-Mayer  potential  correspond  to  the  A  and  b  constants. 
Thus  one  can  tailor  the  frequency  and  hardness  of  collisions 
by  appropriate  selection  of  the  parameters.   At  large  separa- 
tions the  Moliere  and  Born-Mayer  functions  have  similar  shapes 
In  all  of  the  potentials  the  repulsive  wall  dominates  the 
dynamics  of  collisions,  while  the  attractive  well  determines 
whether  an  atom  will  escape  the  surface  and  be  sputtered. 
Reference  [42]  contains  a  detailed  discussion  on  potential 
function  selection. 

2.   Solid  Phase  Potentials 

The  Ti-Ti  interactions  within  the  crystal  were 

o 

modeled  using  a  standard  Moliere-Morse  spline  (a  =  0.10529  A) 
to  allow  for  the  anticipated  hard  collisions  in  the  experi- 

o 

ment.   R  was  set  at  the  lattice  parameter  value  of  2.95  A 
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ad  Dg  was  adjusted  to  give  the  correct  crystal  cohesive 
eergy  of  4.855  eV.   The  Morse  alpha  was  set  at  1.33  A*1 
s>  that  the  slope  closely  matched  that  of  the  Moliere.   The 
tro  were  then  connected  by  a  spline  between  R  =  0.68  and 
I}  =  0.72  LU.   Rc/  the  distance  at  which  the  potential  is 
truncated,  was  set  at  2.20  LU  to  allow  only  interactions  be- 
ween  nearest  and  next-nearest  neighbor  atoms.   Since  the 
orces  produced  at  greater  separations  are  very  weak,  this 
election  of  R  allows  calculation  time  to  be  minimized  with 
egligible  effect  on  the  results. 

The  Nb-Nb  and  V-V  potentials  were  similarly  chosen 

°-l  °-l 

nth  alpha  equal  to  1.34A   and  1.62  A   respectively.   R_ 

/as  set  at  2.2  LU  in  both  cases.   R  was  calculated  as  /3/2 

:imes  the  lattice  constant  for  the  bcc  structure. 

The  metal  (M) -Ar  potential  functions  used  for  the 

bombardment  of  the  crystal  by  argon  were  chosen  as  standard 

Moliere  potentials  since  no  binding  is  to  be  expected.   Rc 

was  set  at  1.7  LU  for  all  runs. 

Morse  potentials  were  used  for  the  M-0  surface 

interactions.   Lacking  any  specific  data  in  the  literature, 

adatoms  in  the  A-top  position  were  placed  at  the  M-0  gas 

phase  equilibrium  separation  distance  above  the  substrate, 

although  this  is  known  to  be  too  high.   The  well  depth  and 

alpha  were  then  adjusted  so  that  the  correct  surface  binding 

energy  was  obtained.   For  the  titanium  two-fold  bridge,  gas 

phase  equilibrium  separation  was  maintained  between  the  oxygen 
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and  its  two  neighbor  titaniums  which  necessitated  placing  the 
oxygens  0.322  LU  above  the  surface  plane.   The  well  depth 
was  again  adjusted  for  the  correct  SBE.   Adatoms  in  the 
three-fold  sites  were  placed  in  the  center  of  the  titanium 
triangles  on  the  surface  plane,  which  required  increasing 

o 

r   to  1.70  A.   R-,  for  all  was  set  just  outside  the  nearest 
e  C  J 

neighbor  distance  to  prevent  0-0  interference.  Parameter 
values  for  the  various  adatom  positions  are  summarized  in 
Table  4 . 

The  Ar-0  potential  used  is  a  Born-Mayer  with  A  =  14.67 

°-l 
keV  and  b  =  -4.593  A   .   This  choice  gave  good  results  in 

previous  work  by  Harrison,  et  al .  [Ref.  21] . 

The  0-0  potential  is  another  Morse  with  DQ  =  5.154  5 

°-l 
eV,  r   =  1.208  A  and  alpha  =  2.691  A   ;  values  which  were 
e 

obtained  from  gas  phase  data  and  were  also  used  by  Harrison, 
et  al . 


TABLE  4 
.  Solid  Phase  Potential  Function  Parameters 


D  (eV) 
e 

re(Al 

Alpha 

(A  1) 

Rc  (1) 

Ti 

ATOP 

5.90 

1.62 

2.0 

2.086 

Ti 

BRIDGE 

2.95 

1.62 

2.0 

2.086 

Ti 

TF 

1.97 

1.70 

2.0 

2.086 

V 

ATOP 

5.52 

1.56 

2.0 

2.584 

Nb 

ATOP 

6.67 

1.69 

2.0 

2.810 
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3.   Gas  Phase  Potentials 

The  potentials  used  for  interactions  above  the 
surface,  as  well  as  in  determining  multimer  formation,  are 
all  of  the  pure  Morse  type.   Parameters  were  selected  from 
References  27,  32  and  33  and  are  summarized  in  Table  5. 
Appendix  A  (Figs.  15-22) ,  contain  plots  of  potential  energy 
versus  separation  for  most  of  the  potential  functions. 

TABLE  5 
Gas  Phase  Potential  Function  Parameters 


-1 


De    (eV) 

re    (A) 

Alpha    (A 

Ti-Ti 

1.466 

2.896 

1.330 

Ti-0 

6.870 

1.620 

2.000 

Nb-Nb 

2.830 

2.858 

1.340 

Nb-0 

7.805 

1.690 

2.000 

V-V 

2.511 

2.622 

1.620 

V-0 

6.462 

1.560 

2.000 

0-0 

5.155 

1.208 

2.691 

F.   MICROCRYSTAL 

Microcrystals  composed  of  676  atoms  were  used  for  each 
computer  run.   The  atoms  are  placed  by  a  lattice  generator 
into  the  proper  hep  or  bec  arrangement,  resulting  in  a 
crystal  15  planes  wide  by  15  planes  long  by  six  planes  deep 
(15x6x15) .   The  first  atom  in  the  crystal  is  number  two, 
one  being  the  ion,  and  is  located  at  (0,0,0,0)  .   Atoms  2-114 
comprise  the  top  layer,  115-226  layer  two,  227-339  layer  three, 
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340-451  layer  four,  452-564  layer  five  and  565-676  layer 
six.   Atom  numbers  above  676  correspond  to  oxygen  adatoms . 

Crystal  size  selection  is  based  on  the  problem  of  confine- 
ment which  is  the  attempt  to  encompass  the  entire  surface 
sensitive  portion  of  the  collision  cascade  within  the  target 
lattice.   As  the  target  size  is  increased  the  yield  will 
asymptotically  approach  a  maximum  value.   Past  this  point 
further  increases  in  size  will  only  affect  computer  run  times. 
The  15x6x15  crystal  size  was  selected  after  making  a  series 
of  runs  on  clean  titanium  using  a  high  yield  impact  point. 
Yield  versus  ion  energy  and  crystal  size  is  depicted  in 
Figures  12  and  13.   As  can  be  seen  containment  is  effectively 
achieved  for  ion  energies  up  to  3.0  keV. 

G.   IMPACT  AREA 

The  impact  area  for  each  simulation  is  a  rectangular 
array  of  104  uniformly  distributed  impact  points  adapted  from 
previous  work  [Ref.  46 J .   Figure  14  is  a  graphical  represen- 
tation of  the  points  and  their  coordinates.   The  actual  length 
and  height  of  the  grid  is  automatically  scaled  by  the  pro- 
gram so  that  its  dimensions  match  those  of  the  interatomic 
spacing  in  the  X  and  Z  directions.   Thus  a  grid  positioned 
on  planes  RBX  =  6.0  and  RBZ  =  6.0  would  have  its  four  corners 
located  on  the  centers  of  the  atoms  at  (6,6),  (8,6),  (6,8) 
and  (8,8) . 

The  impact  points  are  slightly  offset  in  the  X  and  Z 
directions  from  true  center  to  prevent  repetitive  sampling 
of  symmetrically  equivalent  points. 
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Fig.  12.  Yield  versus  crystal  size  at  various  ion  energies. 
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Fig.  13  Yield  versus  crystal  depth. 
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Fig.   14     Standard  impact  zone  with  104  impact  points 
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Independent  data  sets  can  be  generated  by  simply  moving 
the  origin  of  the  grid.   Selection  of  grid  position  is  based 
on  the  surface  symmetry  resulting  from  adatom  density  and 
placement.   Generally  two  sets  of  impact  points  are  re- 
quired to  sample  a  completely  representative  area  of  the 
target  surface.   Surfaces  which  could  be  mapped  into  one 
grid  were  still  run  with  two  sets,  shifted  by  one  plane, 
so  that  all  runs  consisted  of  20  8  impact  points. 

H.   ION  ENERGY 

All  simulations  were  conducted  with  an  ion  energy  equal 
to  1.0  keV.   This  energy  produced  a  reasonable  sputtering 
yield  and  allowed  the  runs  to  be  conducted  on  a  moderate 
sized  crystal.   Additional  runs  were  conducted  on  the 
titanium  C(2x2)TFB  surface  at  0.5  and  2.0  keV,  for  compari- 
son.  This  coverage  produced  the  highest  yield  at  1.0  keV 
and  therefore  the  differences  in  results  as  the  ion  energy 
is  varied  should  be  more  pronounced. 

I.   PROGRAM  MODIFICATIONS 
1.   QRAD 

During  the  early  stages  of  this  investigation 
Harrison  completed  extensive  modifications  to  QRAD  resulting 
in  improved  energy  conservation,  timestep  logic  and  output 
format.   Considerable  effort  was  expended  proofing  and  de- 
bugging the  new  version.   Additional  code  was  also  inserted 
to  identify  registered  pairs,  metal-oxygen  pairs  on  the 
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undamaged  surface,  for  later  use  in  tracking  multimer 
origins. 

2.  AN83 

This  is  the  analyzer  program  used  to  determine  vari- 
ous yield  distributions  and  multimer  data.   Changes  were 
made  to  make  the  program  compatible  with  the  new  output 
format  of  QRAD. 

3.  PLQTAN 

This  program  used  the  output  from  QRAD  to  generate 
graphic  displays  of  the  results.   The  program  was  modified 
to  permit  separate  analyses  of  each  atom  type.   It  was  also 
expanded  to  accommodate  additional  crystal  types.   The  logic 
was  modified  to  allow  determination  of  the  angular  distributions 
from  data  generated  by  the  non-symmetric  sets  of  impact 
points  in  use.   Plot  scaling  routines  were  also  changed  to 
handle  the  additional  crystal  types. 

4.  ADPLOT 

ADPLOT  is  used  to  generate  the  pictorial  representa- 
tion of  adatoms  and  steps  on  the  crystal  surface.   Modifica- 
tions were  made  to  permit  real  time  changes  to  the  scaling 
of  the  plot  size  and  to  adjust  adatom  radii. 
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IV.   RESULTS 

A .   GENERAL 

Detailed  information  on  each  run  has  been  assembled  into 
appendices  following  the  main  body  of  this  thesis.   Included 
are  graphical  displays  of  the  ejection  time  distributions, 
atoms  ejected  per  single  ion  (ASI) ,  energy  distributions, 
spot  patterns,  and  the  ejection  time  differences  between 
multimer  components.   Where  applicable  these  have  been  broken 
down  into  separate  substrate  and  adatom  plots.   Plot  headings 
are  generally  self-explanatory,  indicating  ion  energy,  metal 
type  and  crystal  face  used  for  the  particular  run.   For 
example  "1.0  keV  Ti(001)/Ar  +  Ox  P(2x2)"  indicates  that  the 
data  pertains  to  bombardment  of  the  0  01  face  of  titanium 
by  1.0  keV  argon  ions.   Ox  PC2x2)  denotes  the  presence  of 
oxygen  adatoms  in  a  PC2x2)  configuration.   CM)  or  COX)  are 
used  to  identify  whether  the  plot  is  for  the  substrate  or 
for  the  adatoms  in  the  particular  configuration.   The  plots 
contain  information  based  on  a  standard  run  of  20  8  impact 
points.   On  occasion  data  from  a  multiple  set  of  520  impact 
points  is  included  to  amplify  significant  characteristics. 
These  are  identified  by  the  notation  (IPX51 .   Combined  plots 
of  adatoms  in  the  TFA  and  TFB  positions  are  identified  by 
TF CA+B) .   Discussions  of  the  more  noteworthy  results  of  the 
experiments  are  contained  in  the  remainder  of  this  chapter. 
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B.   YIELD  AND  AS I 

The  yields  from  the  runs  conducted  on  the  clean  crystal 
surfaces  are  summarized  in  Tables  6  and  7.   All  runs  consisted 
of  the  standard  set  of  20  8  impact  points  with  the  exception 
of  the  expanded  run  on  titanium  which  was  composed  of  520 
impact  points.   This  larger  run  was  used  to  test  the  sta- 
tistical validity  of  findings  from  the  smaller  sets.   Of 
particular  interest  is  the  fact  that  the  yield  from  the 
representative  sample  of  208  impact  points  for  titanium 
is  within  1.6%  of  the  yield  from  a  much  larger  run  of  520 
impact  points,  demonstrating  the  validity  of  this  rather  small 
set.   The  significantly  higher  yield  of  argon  from  the  bom- 
bardment of  niobium  and  vanadium  is  to  be  expected  in  view 
of  the  lower  mass  ratio,  0.429  (.versus  0.834  and  0.784  for 
Ti  and  V  respectively)  and  the  increase  in  the  collision 
cross-section  which  is  a  function  of  the  atomic  number. 

The  increased  yield  of  V  and  Nb,  as  compared  to  Ti,  may 
be  explained  by  two  characteristics: 

1.  The  more  open  nature  of  the  bcc  (110)  face.   While 
the  overall  cohesive  energy  is  greater,  the  surface  layer 
has  fewer  nearest  neighbors  and  therefore  fewer  bonds  must 
be  broken  for  an  atom  to  escape  the  surface. 

2.  The  fact  that  there  are  no  channels  present  in  the 

bcc  (.110)  face  for  bombardment  at  normal  angles  of  incidence. 
Hence  all  the  ion  energy  is  deposited  in  the  top  layers  of 
the  crystal  from  which  the  vast  majority  of  sputtered  atoms 
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COVERAGE 


CLEAN 


TABLE  6 
Yield  and  AS I  (by  Coverage) 


NM  (%) 


*U  (%) 


N 


AR 


M 


"0 


Ti  (520) 

945 

- 

10 

1.8 

_ 

Ti  (208) 

374 

- 

2 

1.8 

_ 

V 

946 

- 

8 

4.5 

_ 

Nb 

682 

- 

52 

3.3 

- 

2x2) 

Ti  ATOP 

291 

(79) 

77 

(21) 

2 

1.4 

0.4 

V   ATOP 

616 

(83) 

126 

(17) 

6 

3.0 

0.6 

Nb  ATOP 

522 

(82) 

112 

(18) 

37 

2.5 

0.5 

Ti  BR 

331 

(80) 

82 

(20) 

4 

1.6 

0.4 

Ti  TFA 

324 

(90) 

38 

(10) 

2 

1.6 

0.2 

Ti  TFB 

386 

(87) 

58 

(13) 

4 

1.9 

0.3 

2x2) 

Ti  ATOP 

215 

(62) 

133 

(38) 

2 

1.0 

0.6 

V   ATOP 

457 

(70) 

198 

(30) 

6 

2.2 

0.9 

Nb  ATOP 

374 

(70) 

159 

(30) 

41 

1.8 

0.7 

Ti  Br 

271 

(65) 

148 

(35) 

4 

1.3 

0.7 

Ti  TFA 

290 

(77) 

86 

(23) 

4 

1.4 

0.4 

Ti  TFB(0 

.5) 

233 

(69) 

105 

(31) 

7 

1.1 

0.5 

Ti  TFB(1 

.0) 

368 

(75) 

124 

(25) 

6 

1.8 

0.6 

Ti  TFB (2 

.0) 

419 

(78) 

116 

(22) 

2 

2.1 

0.6 

NOTE:   All  ion  energies  are  1.0  keV  unless  otherwise  noted, 

i.e.,  (0.5).   The  designations  (208)  and  (520) 

refer  to  a  single  set  of  impact  points  and  a 

combined  group  of  520  impact  points  respectively. 

N,   is  the  number  of  reflected  argon  ions.   N,„ 
Ar  3  M 


and  N   are  the  actual  counts  of  ejected  atoms, 
and  Y  are  the  computed  average  yields,  atoms/ion 


'M 
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TABLE  7 
Yield  and  ASI  (by  Substrate) 


COVERAGE 


CLEAN 

Ti 
Ti 


(520) 
(208) 


P(2X2) 

Ti  ATOP 

Ti  BR 

Ti  TF  <A+B> 

CC2X2) 

Ti  ATOP 

Ti  BR 

Ti  TF  <A+B> 


V  CLEAN 

V  P(2X2) 

V  C(2X2) 

Nb  CLEAN 
Nb  P(2X2) 
Nb  C(2X2) 


ATOP 
ATOP 


ATOP 
ATOP 


NM  (%) 


945 
374 


291  (79) 
331  (80) 
355  (88) 


215  (62) 

271  (65) 

329  (76) 

946 

616  (83) 

457  (70) 

682 

522  (82) 

374  (70) 


N 


O 


(%) 


N 


AR 


77  (21) 
82  (20) 
48  (12) 


133  (38) 
148  (35) 
105  (24) 


126 
198 


112 
159 


(17) 
(30) 


(18) 
C30) 


10 
2 


2 
4 
3 


2 

4 
5 

8 
6 
6 

52 
37 

41 


'M 


1.8 
1.8 


1.4 
1.6 
1.8 


1.0 
1.3 
1.6 

4.5 
3.0 
2.2 

3.3 
2.5 
1.8 


0.4 
0.4 
0.3 


NOTE:   The  quantity  <A+B>  for  the  three-fold  sites 

refers  to  the  combined  average  of  the  A  and  B 
positions  of  the  adatoms. 
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originate.   The  yield  per  layer  data  for  adatoms  in  the  A- 

top  position  is  given  in  Table  8.   Yield  ratios  for  the 

other  adatoms  locations  were  found  to  be  essentially  identical 

TABLE  8 
Sputtering  Yield  by  Layer 


COVERAGE 

CLEAN 

Ti  (001) 

V  (110) 
NB  (110) 

P(2x2) 

Ti  (001) 

V  (110) 
Nb  (110) 

C(2x2) 

Ti  (001) 

V  (110) 
Nb  (110) 


LAYER  1 

372 

900 

694 

285 

566 

490 

213 

429 

34  7 

LAYER  2 


0 
44 
26 


5 
41 
24 


2 
24 
24 


LAYER  3 


0 
8 
8 


1 
7 
8 


0 
8 
3 


LAYER  4 


0 
0 
0 


0 

1 

0 


With  respect  to  the  various  adatom  coverages  the  metal 
yields  for  the  A-top  positions  are  fairly  consistent  at 
about  80%  for  the  P(2x2)  and  70%  for  the  C(2x2).   As  the 
oxygen  coverage  is  increased  there  are  fewer  direct  impacts 
on  the  substrate  atoms  which  begin  to  look  more  like  a  second 
layer  to  the  incoming  ions  with  a  corresponding  decrease  in 
yield.   This  decrease  is  summarized  below  where  Delta  Y.. 
represents  the  drop  in  substrate  yield  as  the  coverage  is 
increased  from  a  quarter  to  a  half  monolayer. 
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TABLE  9 
Change  in  Metal  Yield  as  a  Function  of  Adatom  Coverage 


COVERAGE  DELTA  Y   (%) 


Ti  P 

-*■ 

C  ATOP 

-17 

V   P 

->- 

C  ATOP 

-13 

Nb  P 

-> 

C  ATOP 

-12 

Ti  P 

-»■ 

C  BR 

-15 

Ti  P 

-> 

C  TFA 

-13 

Ti  P 

->■ 

C  TFB 

-12 

As  the  adatom  position  is  shifted  from  A-top  to  bridge 
to  three-fold,  the  total  yield  (M+0)  fluctuates  with  the 
maximum  coming  from  the  bridge  positions.   The  metal  yields 
are  about  the  same  percentage  for  both  the  A-top  and  bridge, 
but  increase  in  the  three-fold  positions.   The  corresponding 
decrease  in  the  fraction  of  oxygens  ejected  is  reasonable 
in  light  of  the  added  stability  of  the  triply  bonded  positions. 

The  increased  metal  yield  in  the  bridge,  and  particularly 
in  the  TFB  positions,  can  again  be  explained  in  terms  of 
channel  blockage.   The  bridge  locations  partially  obstruct 
the  channel  in  the  hep  CO 01)  plane,  while  the  TFB  position 
above  the  channel  closes  it  to  incoming  ions. 

The  Atoms  per  Single  Ion  graphs  for  selected  cases  are 
contained  in  Appendix  B.   The  distributions  for  the  substrate 
atoms  are  in  general  broader  than  for  the  adatoms.   As  the 
oxygen  density  is  increased,  there  is  little  change  in  the 
substrate  ASI  distribution,  but  the  maximum  in  the  oxygen 
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distributions  shifts  toward  an  ASI  of  two.   This  is  to  be 
expected  due  to  the  larger  number  of  oxygen  atoms  on  the 
surface.   The  ASI  distributions  of  the  substrate  broaden 
somewhat  as  the  adatom  positions  are  changed  from  the  bridge 
to  the  three-fold,  consistent  with  the  increasing  metal 
yield. 

C.   EJECTED  PARTICLE  ENERGY  VERSUS  COVERAGE 

It  is  interesting  to  note  that  as  the  adatom  coverage 
in  the  A- top  position  is  increased,  the  average  energy  of 
the  ejected  substrate  atoms  also  increases.   As  noted  above, 
this  corresponds  to  a  decrease  in  metal  atom  yield  as  well. 
Lower  energy  particles  are  apparently  retained  in  the  crystal, 
thus  raising  the  average  ejection  energy.   The  most  probable 
energy  of  the  ejected  metal  atoms  also  is  significantly 
higher  when  there  are  adatoms  in  the  A-top  position  and  in- 
creases with  oxygen  density  as  well  (Table  10) .   This  trend 
does  not  extend  to  the  bridge  and  three-fold  positions  in 
titanium,  but  the  most  probable  energies  of  the  substrate 
atoms  are  more  in  line  with  those  ejected  from  clean  metal 
surfaces. 

There  is  a  much  closer  correlation  between  the  metal  and 
oxygen  average  energies  in  the  bridge  and  three-fold  posi- 
tions.  The  multiple  bonding  involved  appears  to  equalize 
the  energies  of  the  ejected  particles  to  some  degree. 

As  the  ion  energy  is  increased  there  is  a  corresponding 
rise  in  the  average  ejection  energy  of  the  metal  atoms,  as 
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TABLE  10 
Average  and  Most  Probable  Ejection  Energies 


COVERAGE 

<EM>  (eV) 

[EMJ  MP  (eV) 

<EQ>  (eV) 

CLEAN 

Ti 
V 

Nb 

10.3 

8.9 

12.3 

3.5 
4.0 

4.5 

~ 

P(2x2) 

(ATOP) 

Ti 
V 

Nb 

10.7 
10.7 
13.8 

7.0 
7.5 
9.0 

11.9 
11.4 
11.9 

C(2x2) 

(ATOP) 

Ti 
V 

Nb 

11.3 
12.2 
15.9 

7.5 
8.0 
9.0 

27.5 
10.7 
13.8 

TITANIUM 

P(2x2) 
C(2x2) 

BR 
BR 

10.057 
9.815 

4.0 
5.0 

10.062 
9.461 

P(2x2) 
C(2x2) 

TFA 
TFA 

9.852 
10.193 

5.0 

4.5 

10.997 
10.190 

P(2x2) 
C(2x2) 

TFB 
TFB 

10.230 
9.743 

3.0 
3.0 

8.055 
9.384 
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TABLE  10 
Average  and  Most  Probable  Ejection  Energies 


COVERAGE 


CLEAN 

Ti 
V 

Nb 

P(2x2)  (ATOP) 

Ti 
V 

Nb 

C(2x2)  (ATOP) 

Ti 
V 

Nb 

TITANIUM 

P(2x2)  BR 
C(2x2)  BR 

P(2x2)  TFA 
C(2x2)  TFA 

P(2x2)  TFB 
C(2x2)  TFB 


<EM>  (eV) 

[EM]  MP  (eV) 

<EQ>  (eV) 

10.3 

3.5 

8.9 

4.0 

— 

12.3 

4.5 

— 

10.7 

7.0 

11.9 

10.7 

7.5 

11.4 

13.8 

9.0 

11.9 

11.3 

7.5 

27.5 

12.2 

8.0 

10.7 

15.9 

9.0 

13.8 

10.057 

4.0 

10.062 

9.815 

5.0 

9.461 

9.852 

5.0 

10.997 

10.193 

4.5 

10.190 

10.230 

3.0 

8.055 

9.743 

3.0 

9.384 
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TABLE  11 
Average  Ejection  Energy  Versus  Ion  Energy 

ION  ENERGY  (keV)         Ti  (eV)  Ox  (eV) 

0.5  8.059  22.514 

1.0  9.743  9.348 

2.0  13.146  20.316 

one  would  expect.   The  changes  in  the  oxygen  energy  match 
the  variation  in  yield  (105,  124,  116) .   While  the  ratio  of 
energies  between  Ti  and  0  is  essentially  the  same  for  0.5 
and  2.0  keV,  the  oxygen  energy  is  about  double  that  for  the 
1.0  keV  ions  at  all  coverages  except  the  A- top  C(2x2)  where 
there  is  also  an  unexpected  rise. 

The  energy  distribution  curves  are  contained  in  Appendix 

C.  The  curves  for  the  ejected  oxygen  atoms  are  typically 
flat  for  the  bridge  and  the  three-fold  positions,  showing 
no  dominant  ejection  energy.   However  there  is  a  pronounced 
Rayleigh-type  shape  to  the  curves  for  oxygen  in  the  A-top 
position  and  they  bear  more  of  a  resemblance  to  those  of  the 
substrate.   This  agrees  with  what  was  found  earlier  where 
the  multiple  bonds  tend  to  cause  a  more  even  distribution  of 
ejected  particle  energy. 

D.  EJECTION  TIME  DISTRIBUTIONS 

Appendix  D  contains  bar  charts  of  the  ejection  time  dis- 
tributions for  representative  runs.   The  distributions  are 
similar  for  titanium  and  vanadium,  but  the  average  ejection 
time  for  niobium  increases  about  30  fsec  to  about  110  fsec 
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(1  femptosecond  =  10~    sec) .   The  distribution  is  also 
broader.   There  is  essentially  no  change  in  any  of  the  dis- 
tributions as  oxygen  is  added  to  the  surface,  or  increased 
from  P(2x2)  to  C(2x2).   Thus  the  ejection  time  appears  to 
be  mainly  a  function  of  the  substrate  mass  and  independent 
of  adatom  coverage  and  lattice  type. 

The  oxygen  ejection  time  distributions  show  the  same 
characteristics  as  the  substrate  atoms.   The  only  significant 
difference  was  that  oxygen  ejection  started  approximately 
10  fsec  later  for  titanium  than  for  vanadium  or  niobium. 
Again  the  distribution  was  broader  for  niobium,  reflecting 
the  influence  of  the  substrate. 

E.   YIELD  VERSUS  ENERGY 

Since  the  C(2x2)  TFB  coverage  produced  the  highest  titanium 
sputtering  yield,  it  was  chosen  to  investigate  the  effects 
of  changing  ion  energy.   This  coverage  was  therefore  run 
at  energies  of  0.5  and  2.0  keV  in  addition  to  the  1.0  keV 
trial.   The  yield  results  and  layer  data  are  contained  in 
Tables  12  and  13. 

TABLE  12 

Titanium  Atoms  Ejected  and  Yield 
as  a  Function  of  Ion  Energy 


ENERGY 
(keV) 

NM 

YM 

No 

Yo 

N7V 

Ar 

0.5 

233 

1.1 

105 

0.5 

0 

1.0 

368 

1.8 

124 

0.6 

0 

2.0 

419 

2.0 

116 

0.6 

0 
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TABLE  13 

Titanium  Atom  Yield  per  Layer 
as  a  Function  of  Ion  Energy 


ENERGY 
(keV) 

LAYER 

1 

LAYER 

2 

LAYER  3 

LAYER  4 

0.5 

232 

1 

0 

0 

1.0 

363 

4 

1 

0 

2.0 

39  8 

16 

1 

4 

As  can  be  seen  the  higher  the  energy,  the  higher  the  yield, 
but  the  relationship  is  definitely  non-linear.   With  increasing 
energy  the  average  ion  penetrates  further  into  the  crystal. 
As  a  consequence  a  smaller  percentage  of  its  energy  is  deposited 
in  the  top  layers.   More  energy  is  carried  deep  into  the 
crystal  where  it  is  dissipated  through  collision  cascades 
producing  only  lattice  displacements  and  vibrations. 

The  decrease  in  oxygen  yield  for  the  2.0  keV  ions  is 
also  probably  caused  by  the  deeper  deposition  of  energy. 
The  energy  distribution  plot  for  the  ejected  oxygens  (Figs. 
54-56)  clearly  shows  a  decrease  in  the  energy  as  compared  to 
the  1.0  keV  case.   Hence  fewer  adatoms  are  receiving  enough 
energy  to  surmount  the  surface  binding  well  and  escape. 

F.   MULTIMER  YIELD 

Tables  14  and  15  present  the  various  multimer  yields  in 
terms  of  total  numbers  and  as  percentages.   Tables  16  and  17 
contain  a  breakdown  by  species  arranged  by  substrate  type 
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TABLE  14 
Multimer  Yield  as  a  Function  of  Coverage 


COVERAGE 


Clean 


1-MER    2-MER    3-MER    4-MER 


Ti  (520) 

860 

43 

Ti  (208) 

356 

10 

V 

731 

74 

Nb 

585 

48 

2x2) 

Ti  ATOP 

304 

25 

V  ATOP 

529 

77 

Nb  ATOP 

441 

60 

Ti  BR 

350 

26 

Ti  TFA 

337 

10 

Ti  TFB 

389 

23 

2x2) 

Ti  ATOP 

250 

49 

V  ATOP 

416 

99 

Nb  ATOP 

377 

59 

Ti  BR 

335 

31 

Ti  TFA 

336 

16 

Ti  TFB  (  .5) 

310 

14 

Ti  TFB  (1.) 

437 

18 

Ti  TFB  (2.) 

490 

17 

2 

2 

22 

13 


4 
12 
18 

3 

1 
1 


2 

3 
7 
3 
2 
0 
5 
3 


0 
0 
2 
2 


0 
3 
2 
2 
0 
0 
0 
0 


TOTAL 

MER 

MULTIMER 

0 

45- 

0 

12 

0 

98 

0 

63 

0 

29 

0 

94 

0 

80 

0 

29 

0 

11 

0 

26 

0 

51 

2 

107 

0 

68 

0 

36 

0 

18 

0 

14 

0 

23 

0 

20 

NOTE:   All  ion  energies  are  1.0  keV  except  as  indicated, 
i.e. ,  (0.5)  . 
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TABLE  15 
Multimer  Yield  (Percentage) 


%  YIELD 

THAT  FORMED 

COVERAGE 

1-MER 

2-MER 

3-MER 

4-MER 

5-MER 

MULTIMERS 

Clean 

Ti  (520) 

90.4 

9.0 

0.6 

0.0 

0.0 

9.6 

Ti  (208) 

97.9 

5.3 

1.6 

0.0 

0.0 

6.9 

V 

76.9 

15.4 

6.9 

0.8 

0.0 

23.1 

Nb 

80.5 

13.1 

5.3 

1.1 

0.0 

19.5 

P(2x2) 

Ti  ATOP 

83.2 

13.5 

3.3 

0.0 

0.0 

16.8 

V  ATOP 

71.6 

20.8 

4.9 

2.7 

0.0 

28.4 

Nb  ATOP 

71.3 

18.9 

8.5 

1.3 

0.0 

28.7 

Ti  BR 

85.2 

12.6 

2.2 

0.0 

0.0 

14.8 

Ti  TFA 

93.7 

5.5 

0.8 

0.0 

0.0 

6.3 

Ti  TFB 

87.1 

10.4 

0.7 

1.8 

0.0 

12.9 

C(2x2) 

Ti  ATOP 

70.1 

28.2 

1.7 

0.0 

0.0 

29.9 

V  ATOP 

65.1 

30.2 

1.4 

1.8 

1.5 

34.9 

Nb  ATOP 

72.5 

22.1 

3.9 

1.5 

0.0 

27.5 

Ti  BR 

82.2 

14.8 

2.1 

1.9 

0.0 

18.8 

Ti  TFA 

89.9 

8.5 

1.6 

0.0 

0.0 

10.1 

Ti  TFB  ( .5) 

91.7 

8.3 

0.0 

0.0 

0.0 

8.3 

Ti  TRB  (1.) 

89.6 

7.3 

3.1 

0.0 

0.0 

10.4 

Ti  TFB  (2.) 

93.6 

6.4 

1.7 

0.0 

0.0 

8.1 

NOTE:   All  ion  energies  are  1.0  keV  except  as  indicated, 
i.e. ,  (0.5). 
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TABLE  16 
Multimer  Species  Formed  (by  Substrate) 


COVERAGE  M2    M3    M4    MO    M20    M30    M02    MO         M302 

Clean 

Ti  (520)  43     20- 

Ti  (208)  10     20- 

P(2x2) 

Ti  ATOP  300     22     4      0      0      0      0 

Ti  BR  8     0    0     18     3      0      0      0      0 

Ti  TF  <A+B>  700      81      1      0      0      0 

C(2x2) 

Ti  ATOP  000492      0      0      0      0 

Ti  BR  5     0    0     26     2      1      1      1      0 

Ti  TF  <A+B>  610     11     2      0      1      0      0 

V  CLEAN  74    22    2     - 

V  P(2x2)  18     0    2     59    12      3      0      0      0 

V  C(2x2)  900     90     3      1      0      2      2 

Nb  CLEAN  48    13    2- 

Nb  P(2x2)  18     6    0     42    12      1      0      1      0 

NB  C(2x2)  220     57     5      0      1      2      0 

NOTE:   The  quantity  <A+B>  refers  to  the  average  of  the  two 
three- fold  positions. 
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TABLE  17 
Multimer  Species  Formed  (by  Coverage) 


COVERAGE 


ML 


M. 


M, 


xMO    M20   MO 


MO. 


M2°2 


M3°2 


Clean 


Ti  (520) 

43 

2 

0 

— 

— 

_ 

_ 

_ 

— 

Ti  (208) 

10 

2 

0 

— 

_ 

_ 

_ 

_ 

_ 

V 

74 

22 

2 

— 

_ 

_ 

_ 

_ 

_ 

Nb 

48 

13 

2 

- 

- 

- 

- 

- 

- 

2x2) 

Ti  ATOP 

3 

0 

0 

22 

4 

0 

0 

0 

0 

V  ATOP 

18 

0 

2 

59 

12 

3 

0 

0 

0 

Nb  ATOP 

18 

6 

0 

42 

12 

1 

0 

1 

0 

Ti  BR 

8 

0 

0 

18 

3 

0 

0 

0 

0 

Ti  BFA 

6 

0 

0 

4 

1 

0 

0 

0 

0 

Ti  TFB 

11 

0 

0 

12 

1 

2 

0 

0 

0 

2x2) 

Ti  ATOP 

0 

0 

0 

49 

2 

0 

0 

0 

0 

V  ATOP 

9 

0 

0 

90 

3 

1 

0 

2 

2 

Nb  ATOP 

2 

2 

0 

57 

5 

0 

1 

2 

0 

Ti  BR 

5 

0 

0 

26 

2 

1 

1 

1 

0 

Ti  TFA 

6 

0 

0 

10 

2 

0 

0 

0 

0 

Ti  TFB  ( .5) 

4 

0 

0 

10 

0 

0 

0 

0 

0 

Ti  TFB  (1.) 

7 

1 

0 

11 

2 

0 

2 

0 

0 

Ti  TFB  (2.) 

9 

0 

0 

8 

2 

0 

1 

0 

0 

NOTE:   All  ion  energies  are  1.0  keV  except  as  indicated, 
i.e. ,  (0.5) . 
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and  adatom  coverage.   The  multimer  yield  is  seen  to  increase 
with  oxygen  coverage  and  is  due  primarily  to  the  ejection 
of  metal  atoms  with  their  attached  adatoms .   There  is  a  good 
correspondence  in  the  multimer  yield  between  vanadium  and 
niobium  in  the  clean  and  P(2x2)  coverage  which  then  degener- 
ates as  the  coverage  is  increased  to  C(2x2) .   This  can  be 
attributed  to  the  decreasing  influence  of  the  crystal  lattice 
and  increasing  effect  of  the  oxygen  on  the  ejection  mechan- 
isms.  This  will  be  covered  in  greater  detail  in  the 
following  section. 

G.   MULTIMER  ORIGINS 

One  of  the  primary  objectives  was  to  study  the  origins 
of  multimer  components  to  compare  fragmentation  and  recom- 
bination as  modes  of  molecule  formation.   Table  18  presents 
data  as  to  whether  the  constituent  atoms  were  nearest-neighors 
in  the  crystal  prior  to  ejection,  or  in  the  case  of  trimers 
or  greater,  whether  a  portion  of  the  molecule  was  joined. 
In  the  table  the  letter  'C'  represents  the  number  of  multimers 
whose  components  were  coupled  together  on  the  surface;  p — 
partially  connected  multimers,  N — those  multimers  whose 
components  were  not  joined  before  ejection,  %1 — the  percentage 
that  was  formed  by  fragmentation  and  %2 — the  percentage 
that  were  formed  either  through  total  or  partial  fragmentation, 
where  a  portion  of  the  multimer  was  joined  on  the  surface. 
Table  19  contains  similar  data  for  trimers  and  higher  order 
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TABLE  18 

Fragmentation  versus  Recombination 
(All  multimers) 


COVERAGE 


Clean 


N 


%1 


%2 


Ti  (520) 

29 

Ti  (208) 

6 

V 

71 

Nb 

41 

2x2) 

Ti  ATOP 

24 

V   ATOP 

76 

Nb  ATOP 

52 

Ti  BR 

23 

Ti  TFA 

9 

Ti  TFB 

21 

2x2) 

Ti  ATOP 

49 

V  ATOP 

86 

Nb  ATOP 

49 

Ti  BR 

29 

Ti  TFA 

11 

Ti  TRB  (.5) 

10 

Ti  TFB  (1.) 

17 

Ti  TFB  (2.) 

13 

1 

15 

64.4 

66.6 

1 

5 

50.0 

58.3 

7 

20 

72.4 

79.6 

1 

21 

65.1 

66.7 

0 

5 

82.8 

82.8 

6 

12 

80.9 

87.2 

7 

21 

65.0 

73.8 

2 

4 

79.3 

86.2 

1 

1 

81.8 

90.9 

0 

5 

80.8 

80.8 

0 

2 

96.0 

96.0 

2 

19 

80.4 

82.2 

4 

15 

72.1 

77.9 

2 

5 

80.6 

86.1 

0 

7 

61.6 

61.6 

0 

4 

71.4 

71.4 

1 

5 

73.9 

78.3 

1 

6 

65.0 

70.0 

NOTE:   C  represents  the  number  of  multimers  that  were  nearest 
neighbors  on  the  surface  prior  to  ejection;  P,  those 
partially  connected  and  N,  the  number  formed  by 
recombination  of  non-adjacent  atoms.   %1  is  the 
percent  formed  by  fragmentation,  %2  the  percent  formed 
by  partial  fragmentation.   All  ion  energies  are  1.0 
keV  except  as  indicated,  i.e.,  (0.5). 
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TABLE  19 

Fragmentation  versus  Recombination 
(Trimers  and  Greater) 


COVERAGE 


N 


%1 


Clean 


Ti  (520) 

1 

Ti  (208) 

1 

V 

17 

Nb 

14 

2x2) 

Ti  ATOP 

23 

V  ATOP 

11 

Nb  ATOP 

13 

Ti  BR 

1 

Ti  TFA 

0 

Ti  TFB 

3 

2x2) 

Ti  ATOP 

1 

V  ATOP 

6 

Nb  ATOP 

5 

Ti  BR 

4 

Ti  TRA 

1 

Ti  TRB  (1.) 

3 

Ti  TFB  (2.) 

1 

1 
1 

7 
1 


0 
6 
7 
2 
1 
0 


0 
2 

4 
2 
0 

1 
1 


0 

50.0 

100.0 

0 

50.0 

100.0 

0 

70.8 

100.0 

0 

93.3 

100.0 

3 

88.5 

88.5 

0 

64.7 

100.0 

0 

65.0 

100.0 

0 

33.3 

100.0 

0 

0.0 

100.0 

0 

100.0 

100.0 

1 

50.0 

50.0 

0 

75.8 

100.0 

0 

55.6 

100.0 

0 

60.0 

100.0 

1 

50.0 

50.0 

1 

60.0 

80.0 

1 

33.3 

66.7 

NOTE:   C  represents  the  number  of  multimers  that  were 
nearest  neighbors  on  the  surface  prior  to  ejec^ 
tion;  P,  those  partially  connected  and  N,  the 
number  formed  by  recombination  of  non-adjacent 
atoms.   %1  is  the  percent  formed  by  fragmenta- 
tion, %2  the  percent  formed  by  partial  fragmenta- 
tion.  All  ion  energies  are  1.0  keV  except  as 
indicated,  i.e.,  (0.5). 
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multimers.   Table  20  presents  the  percentage  of  multimer 
species  formed  by  fragmentation. 

The  general  trend  is  toward  more  fragmentation  as  the 
oxygen  coverage  is  increased.   This  is  due  primarily  to  the 
fact  that  substrate  atoms  are  ejected  after  being  hit  from 
underneath  and  carry  off  their  attached  adatoms  as  they  leave. 
Many  partial  fragmentation  molecules  were  formed  when  attached 
atoms  exited  the  surface  obliquely,  knocking  off  oxygens  from 
nieghboring  atoms  on  the  way  out.   This  also  appears  to  be 
one  of  the  principal  mechanisms  involved  in  recombination 
of  non-a joining  atoms. 

Figures  74-87  of  Appendix  E  and  Table  21  give  the  loca- 
tions and  frequency  of  occurrence  for  all  of  the  observed 
mechanisms  involved  in  multimer  formation.   (In  the  table 
double  entries  refer  to  the  a  and  b  sections  of  the  particu- 
lar drawing.)   By  far  the  greatest  source  of  multimers  is 
the  dimer  composed  of  a  substrate  atom  and  its  attached 
adatom  which  is  a  fragmentation  mechanism.   Higher  order 
multimers  are  a  fairly  rare  occurrence,  but  over  65%  are 
formed  by  fragmentation  and  virtually  all  involve  at  least 
partial  fragmentation  (Table  20) . 

As  the  oxygen  coverage  is  increased,  the  monoxide  is 
still  the  prevalent  species  formed,  but  there  is  a  shift  in 
the  mechanisms  that  make  up  the  rest  of  the  yield  and  new 
mechanisms  are  added.   In  the  A- top  and  bridge  positions 
the  percentage  of  multimers  resulting  from  fragmentation 
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TABLE  20 

Percentage  of  Multimer  Species 
Formed  by  Fragmentation 


COVERAGE 


M. 


M. 


M, 


MO 


M20 


M30 


MO. 


M  O 
2U2 


M3°2 


Clean 


Ti 

67 

50 

- 

- 

— 

— 

— 

— 

V 

74 

73 

100 

— 

— 

_ 

_ 

_ 

Nb 

56 

92 

100 

- 

- 

- 

- 

- 

2x2) 

Ti  ATOP 

33 

— 

— 

82 

100 

_ 

_ 

_ 

V  ATOP 

100 

- 

100 

80 

50 

100 

- 

- 

Nb  ATOP 

83 

100 

- 

57 

50 

100 

— 

IP 

Ti  BR 

63 

- 

- 

89 

33 

— 

— 

— 

Ti  TFA 

100 

- 

- 

75 

IP 

- 

— 

- 

Ti  TFB 

64 

- 

- 

92 

100 

100 

- 

- 

Ti  TF(A+B) 

82 

- 

- 

84 

50 

100 

- 

- 

2x2) 

Ti  ATOP 

- 

— 

_ 

98 

50 

— 

_ 

— 

V   ATOP 

56 

- 

- 

81 

67 

IP 

- 

100 

Nb  ATOP 

100 

100 

- 

74 

40 

- 

X 

X 

Ti  BR 

40 

- 

- 

92 

100 

100 

IP 

IP 

Ti  TRA 

33 

- 

- 

80 

50 

— 

- 

- 

Ti  TFB 

86 

X 

- 

73 

100 

— 

50 

- 

Ti  TR(A+B) 

60 

X 

- 

77 

75 

— 

50 

- 

100 


/ 

NOTE:   A  tack  (-)  indicates  that  the  species  was  not  formed, 
while  "X"  indicates  that  none  were  formed  by  frag- 
mentation.  IP  indicates  that  only  one  of  the  particu- 
lar species  was  formed  and  that  was  through  partial 
fragmentation . 
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TABLE    21 

Frequency  of  Occurrence  of   the  Various 
Ejection  Mechanisms 


COVERAGE 


8     9     10     11     12 


dean 


Ti  (520) 

20/8 

1 

2 

9/3 

1 

1 

— 

Ti  (208) 

3/2 

2 

2 

1 

- 

- 

— 

— 

— 

— 

_ 

— 

V 

37 

13 

4 

9 

5 

7 

— 

4 

2 

2 

4 

— 

Nb 

19 

8 

10 

4 

6 

2 

8 

- 

2 

- 

- 

2 

[2x2) 

Ti  ATOP 

1 

1 

19 

1/2 

2 

— 

— 

— 

_ 

_ 

_ 

_ 

V  ATOP 

6/0 

6 

1 

- 

- 

2 

47/10 

- 

- 

4 

2 

3 

NB  ATOP 

5/1 

3 

3 

3 

4 

- 

24/7 

7 

1 

5 

1 

2 

Ti  BR 

4 

1 

1 

1 

1 

12 

2 

1 

1 

— 

— 

— 

Ti  TEA 

1 

2 

1 

1 

3 

1 

— 

Ti  TFB 

1 

5 

1 

1 

1 

1 

1 

11 

1 

1 

- 

- 

!2x2) 

Ti  ATOP 

1 

45 

1 

1 

— 

— 

_ 

_ 

_ 

— 

— 

_ 

V  ATOP 

3 

1 

9 

1 

2 

- 

74 

5 

7 

3 

- 

- 

Nb  ATOP 

2 

- 

- 

- 

- 

2 

42 

10 

2 

- 

1 

2 

Ti  BR 

1 

1 

1 

1 

24 

2 

1 

1 

1 

1 

1 

1 

Ti  TRA 

1. 

3 

1 

1 

1/1 

6/2 

1 

1 

- 

- 

- 

- 

Ti  TFB  (.5) 

- 

- 

1/1 

- 

1 

5/2 

- 

- 

- 

- 

- 

1 

Ti  TFB  CI.) 

1 

1/1 

1/0 

2 

1 

4/4 

1/1 

1 

1 

1 

1 

1 

Ti  TFB  (2.) 

1 

- 

3/2 

- 

- 

5/1 

- 

1 

- 

- 

- 

1 

NOTE:     The  column  headings  refer  to  surface  locations  shown  in 

Appendix  E.  Columns  with  double  entries  refer  to  the  a  and 
b  sections  of  the  individual  figures.  Ion  energies  are  1.0 
keV  except  as  indicated,   i.e.,    (.5) . 
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TABLE    21     (CONT.) 


COVERAGE  13     14     15     16     17     18     19     20     21     22     23     24 

Clean 

Ti   (520)  ------------ 

Ti   (208)  ------------ 

V  --2241------ 

Nb  li--______„_ 

P(2x2) 

Ti  ATOP  ____________ 

V  ATOP  21-12-3/2     11--- 
Nb  ATOP  --2--1---123 
TiBR  ____________ 

Ti  TEA  _____________ 

Ti  TFB  ____________ 

C(2x2) 

Ti  ATOP  ____________ 

V  ATOP  -___-___7_i2 
Nb  ATOP  -1121111-1-- 
TiBR  ____________ 

Ti  TEA  ____________ 

Ti  TFB   (.5)        1       ----------- 

Ti  TFB   (1.)        --_-         ________ 

Ti  TFB    (2.)        -321         i---_-_- 


COVERAGE  25     26     27     28 

P(2x2) 

Nb  ATOP  12       11 


76 


increases,  while  it  decreases  for  the  three- fold  positions. 
Those  oxygens,  sharing  three  bonds,  tend  to  be  left  behind 
unless  at  least  two  of  their  nearest  neighbors  are  also 
ejected.   The  greater  oxygen  density  also  presents  a  larger 
target  to  obliquely  ejected  fragments,  increasing  the 
probability  of  collision  and  subsequent  capture. 

H.   MULTIMER  EJECTION  TIMES 

Appendix  F  contains  bar  charts  depicting  the  differences 
in  ejection  times  between  multimer  components,  both  for 
metal-metal  and  metal-oxygen  pairs.   Essentially  none  of 
the  pairs  in  any  of  the  runs  were  ejected  simultaneously, 
but  all  but  about  seven  pairs  were  ejected  within  100  fsec 
of  each  other  and  most  within  50  fsec.   This  difference  is 

about  the  same  order  of  magnitude  as  the  vibrational  period 

-13 
of  an  equivalent  diatomic  molecule  [Ref .  4  7] ,  10     sec 

(10  0  fsec) ,  and  is  not  necessarily  an  indication  that  the 
atoms  were  ever  disjoined  at  some  point  in  time. 

Comparison  of  the  charts  for  the  different  oxygen  cover- 
ages reveals  a  correspondence  between  ejection  time  differ- 
ence and  the  percentage  of  multimers  that  were  formed  by 
fragmentation.   As  one  would  expect  there  is  an  increasing 
time  difference  between  component  ejection  as  the  number  of 
molecules  formed  by  recombination  increases. 

I.   SPOT  PATTERNS 

The  spot  patterns  produced  by  ion  bombardment  of  the 
crystal  surfaces  are  contained  in  Appendix  G.   The  ejection 
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patterns  correspond  to  the  directions  of  the  close-packed 
rows  in  the  crystal,  which  for  hep  are  oriented  at  30,  90 
and  150  degrees  from  the  <001>  direction.   For  bec  the  direc- 
tions are  0  and  90  degrees.   These  lead  to  a  six- fold  symmetry 
pattern  in  hep  and  a  four-fold  rectangular  pattern  in  bcc. 
The  patterns  for  vanadium  and  niobium  are  in  quite  good 
agreement  with  the  expected  results,  but  there  is  consider- 
able deviation  in  the  case  of  titanium  atoms.   The  spikes  at 
90,  150  and  330  degrees  are  present,  but  there  are  extra 
spikes  at  0,  60,  180  and  240  degrees,  while  those  at  30  and 
210  degrees  are  missing.   A  variety  of  different  runs  were 
conducted  with  additional  impact  point  origins  and  permutated 
impact  point  sets,  but  all  resulted  in  the  same  pattern.   At 
present  there  is  no  apparent  reason  for  this  discrepancy. 
Second  layer  contributions  would  tend  to  broaden  the  spikes 
somewhat,  but  not  change  the  overall  appearance. 

The  patterns  for  the  metals  with  adatoms  in  the  A-top 
position  are  essentially  the  same  as  for  those  from  clean 
surfaces,  even  though  they  have  considerably  more  energy  as 
was  pointed  out  earlier.   There  is  some  blurring,  presumably 
due  to  adatom  interference,  but  it  is  very  minor.   This 
blurring  becomes  more  pronounced  as  the  oxygen  density  is 
increased  from  P  to  C.   The  features  of  the  metal  spot 
patterns  begin  to  fade  when  adatoms  are  placed  in  the  bridge 
positions  and  are  completely  absent  for  oxygen  in  the  three- 
fold sites,  leaving  a  uniform  angular  distribution. 
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There  is  no  marked  tendency  for  oxygen  to  eject  perpendicu- 
lar to  the  surface  in  any  of  the  locations,  as  was  found  by 
Hanson  [Ref .  37] .   The  ejection  angles  were  observed  to  in- 
crease slightly  with  increasing  oxygen  density,  probably 
caused  by  the  interference  of  other  adatoms  and  increasing 
yield.   This  was  also  found  to  be  the  case  for  varying  ion 
energy,  where  the  average  ejection  angle  reached  a  maximum 
for  1.0  keV  bombardment  corresponding  to  the  peak  in 
oxygen  yield. 

J.   EJECTION  TIME  VERSUS  ENERGY 

The  plots  of  ejection  time  versus  energy  are  contained  in 
Appendix  H.   As  can  been  seen,  the  highest  energy  atoms  are 
ejected  very  early  in  the  process.   As  time  progresses,  the 
ejection  energy  falls  rapidly  with  a  few  low  energy  stragglers 
coming  off  as  the  collision  cascade  dies  and  recrystalliza- 
tion  begins.   The  distributions  are  independent  of  either 
adatom  density  or  location.   The  duration  of  the  ejection 
period  is  again  seen  to  be  a  function  of  the  crystal  atomic 
mass  as  was  found  earlier  in  comparing  ejection  time  dis- 
tributions.  Increasing  ion  energy  increases  the  energy  of 
those  atoms  which  are  ejected  early,  but  has  no  appreciable 
effect  on  the  bulk  of  the  ejected  atoms. 

The  energy  distribution  of  ejected  oxygens  mimics  that 
of  the  substrate,  with  the  most  energetic  atoms  coming  off 
early  in  the  cascade.   This  close  correlation  is  most  likely 
due  to  fragmentation  where  oxygens  are  ejected  with  their 
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neighboring  metal  atom(s) .   There  is,  though,  a  faster  drop- 
off in  energy  than  in  the  case  of  the  substrate.   Some 
broadening  is  observed  with  more  energetic  atoms  being 
ejected  later  as  the  adsorption  sites  are  shifted  to  those 
involving  a  higher  degree  of  bonding. 

K.   SURFACE  LAYER  EJECTION  FREQUENCY 

Appendix  I  contains  diagrams  depicting  the  relative 
frequency  of  ejection  of  surface  layer  atoms.   For  simplicity, 
the  plots  contain  information  for  a  single  set  of  104  impact 
points.   Also  indicated  are  the  locations  of  the  impact 
zones  that  produced  the  ejections. 

The  most  obvious  feature  is  the  tendency  for  the  damage 
area  to  extend  vertically  in  the  Z  direction  for  the  hep 
crystals,  while  the  pattern  for  bec  is  much  more  symmetric. 
This  can  be  explained  in  terms  of  the  locations  of  close- 
packed  rows.   The  two  close-packed  rows  in  bec  are  in  the  X 
and  Z  directions,  however  hep  only  has  one  of  three  in  the 
X  direction,  with  the  other  two  primarily  in  the  Z  direction. 

The  plots  also  graphically  demonstrate  the  decrease  in 
substrate  yield  as  the  oxygen  density  is  increased,  as  was 
mentioned  in  Section  B. 

L.   COMPARISON  WITH  EXPERIMENTAL  DATA 

Tables  22  and  23  provide  a  comparison  of  experimental 
data  on  the  ratios  of  multimer  species  formed  by  recombina- 
tion and  fragmentation  with  that  obtained  by  computer 
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TABLE    22 

Comparison  of   Experimental  and  Simulation   Results 
on  Multimer  Formation  Mechanisms   for  Titanium 


PCPCPCPC 
RATIO  YU       (2x2)       (2x2)       (2x2)      (2x2)      (2x2)       (2x2)       (2x2)       (2x2) 

ATOP       ATOP         BR  BR         TFA+B     TFA+B       ALL         ALL 


O'/T^         17.5     20.0       50.6       21.0       51.5       11.7       30.3       15.7       38.4 


TiO+/Ti+  50.0   1.55   0.60   0.70   0.92   0.43   1.43   0.70   1.07   R 
TiO"/Ti+   0.01  7.05  29.6   5.40   10.7    1.97   2.16   4.00   8.63   F 


TioJ/Ti*   0.67    -    -     -     0.04     -    0.12 


0.04   R 


TiO~/Ti+   0.28  -     0.36     -    0.18    -     0.26   F 


TiO+/0"    286.  7.80   1.20  3.20    1.80   3.78   4.64   4.47   2.78   R 
TiO"/0"    0.07  35.3   58.5  26.3    20.8   17.2   6.96   25.3   22.5   F 


TiO*/0~    3.90    -         -     0.08         0.44         0.09   R 


'2' 
TiO~/0~    1.60 


0.82     -    0.66    -     0.71   F 


NOTE:   All  quantities  are  in  percentages.   Experimental  data 
taken  from  Reference  22.   "R"  and  "F"  refer  to 
whether  the  mechanism  involved  was  recombination  or 
fragmentation.   The  TFA+B  entries  were  obtained  by 
averaging  the  data  from  the  TFA  and  TFB  positions. 
The  ALL  entires  represent  the  average  from  all  three 
adatom  sites. 
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TABLE    23 


Comparison   of  Experimental   and  Simulation  Results   on 
Multimer  Formation  Mechanisms   for  Vanadium  and  Niobium 


PC  PC 
RATIO             YU            (2x2)           (2x2)                RATIO             YU            (2x2)  (2x2) 
ATOP           ATOP                                                         ATOP           ATOP 


0"/V+         10.0         10.8  29.1  0~/Nb+         UNK         13.6  30.8 


VO+/V+         27.7  2.44  5.19  NbO+/Nb+       UNK  4.55  5.01       R 

VO~A+         0.02  9.76  22.1  NbO~/Nb+       UNK  6.04         14.3         F 


TO^/V+         0.13  -  -  NbO^/Nb"1"       UNK  -  0.42         R 

VO~A+         0.78  -  -  Nbol/Nb*       UNK  -  0.28         F 


VO+/0~         278.  22.6  17.8  NbO+/0~         UNK         32.7  16.6         R 

VO"/0"         0.19  90.4  76.0  NbO~/0~         0.72       43.3  47.4         F 


VO^/O"         0.14  -  -  NbO^/0~         UNK  -  1.32         R 


VO~/0~         7.80  -  -  NbO^/o"         14.0 


0.88         F 


NOTE:      All   quantities   are   in  percentages.      Experimental   data 
taken   from  Reference   23.       "R"    and    "F"    refer  to  whether 
the  mechanism  involved  was  recombination  or 
fragmentation . 
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simulation.   The  experimental  data  was  obtained  by  Yu 
[Refs.  22,23]  using  500  ev  neon  and  argon  to  bombard  the 
three  metals.   His  experiments  used  cesiated-oxygenated 
polycrystals  and  through  the  change  in  the  work  function 
induced  by  the  cesium,  he  was  able  to  obtain  information 
on  the  charge  state  dependence  of  the  two  multimer  ejection 
mechanisms.   The  ion  ratios  presented  in  the  tables  were 
calculated  by  extrapolating  the  relative  yield  of  each 
species  back  to  a  zero  change  in  the  work  function  which 
should  represent  the  cesium-free  case  (Appendix  J) . 

No  exact  comparison  can  be  made  due  to  the  uncertainty 
with  respect  to  relative  ionization  probability,  oxygen 
placement  and  density,  as  well  as  differences  in  ion  mass 
and  energy.   However,  if  the  ratios  of  ions  formed  by  the 
two  mechanisms  were  similar,  the  computer  simulations  would 
support  Yu's  proposed  model.   This  was  not  found  to  be  the 
case.   In  each  of  the  oxygen  positions  and  densities  the 
ratio  of  fragmentation  to  recombination  was  found  to  be 
opposite  that  reported  by  Yu. 

This  variance  between  the  simulation  and  experimental 
results  is  probably  not  due  to  differences  in  ion  energy 
since  Table  12  indicates  that  the  total  multimer  yield 
decreases  as  the  ion  energy  increases.   The  percentage  of 
multimers  formed  by  fragmentation  also  decreases,  but  this 
is  believed  due  to  the  longer-ranged  transport  of  energy 
which  allows  more  ejection  mechanisms  to  operate.   The 
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activation  of  more  mechanisms  could  then  decrease  the  rela- 
tive fraction  of  multimers  formed  by  fragmentation.   The 
simulation  yield  obtained  at  500  eV  was,  in  any  case,  too 
small  for  meaningful  comparison. 

Failure  of  the  simulations  to  support  Yu's  theory  may 
be  caused  by  an  inadequate  treatment  of  the  ionization  proba- 
bility contribution  to  the  determination  of  multimer  sta- 
bility.  Since  the  program  does  not  determine  charge  state, 
some  of  the  multimers  that  were  classified  as  stable  may, 
in  reality,  be  unstable  and  thus  reduce  the  percentage 
found  to  be  formed  by  fragmentation. 
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V.   CONCLUSIONS 

The  program  has  the  inherent  flexiblity  to  permit  de- 
tailed analysis  of  virtually  any  aspect  of  sputtering.   The 
ease  with  which  parameters  may  be  isolated  and  varied 
greatly  facilitates  the  determination  of  the  significant 
factors  involved,  as  well  as  the  interdependence  of  these 
factors. 

The  results  of  these  simulations  support  earlier  find- 
ings concerning  the  effects  of  channeling  and  that  the  majority 
of  sputtered  atoms  originate  in  the  top  layers  of  the  crystal. 
The  angular  distribution  of  ejected  atoms  agrees  with  ex- 
perimental results  in  the  case  of  vanadium  and  niobium, 
however  the  discrepancies  noted  in  the  ejection  patterns  for 
titanium  require  further  study. 

The  decrease  in  oxygen  yield  from  crystals  with  adatoms 
in  the  three-fold  positions  indicates  that  the  oxygens  are 
most  stable  in  these  sites.   Such  positions  would  also 
facilitate  their  incorporation  into  the  metal  lattice  as  the 
oxygen  exposure  is  increased,  resulting  in  the  formation  of 
a  stable  oxide  in  the  top  layers.   It  is  also  possible  that 
the  oxygens  assume  an  A-top  position  at  low  oxygen  exposures 
and  as  the  coverage  is  increased,  oxygen-oxygen  repulsion 
causes  a  migration  to  the  three-fold  sites. 

The  disagreement  between  the  simulation  results  and 
Yu's  theory  and  findings  on  multimer  formation  need  not  be 
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due  to  a  failure  of  his  basic  model,  but  rather  to  an  incom- 
plete understanding  and  treatment  of  the  chemistry  involved 
after  the  atoms  are  ejected  from  the  surface  by  the  simulation 
model.   On  the  other  hand  it  is  somewhat  difficult  to  under- 
stand how  the  yield  of  TiO  ,  presumably  formed  by  recombination 
as  suggested  by  Yu,  can  be  so  much  larger  than  TiO~,  from 
fragmentation,  considering  the  probabilities  involved  in  two 
or  more  atoms  finding  themselves  in  the  same  location  with 
relative  velocities  favorable  enough  to  permit  bonding.   It 
may  also  be  that  a  significant  percentage  of  the  negative 
ions  stabilize  themselves  by  electron  ejection  prior  to 
their  collection  in  the  mass  spectrometer.   Yu  has  also 
pointed  out  that  several  discrepancies  in  his  findings  may 
be  due  to  the  fact  that  as  the  oxygen  density  is  increased 
TiO   and  Ti02  may  be  formed  by  both  mechanisms. 

An  item  that  bears  further  study  is  the  increase  in  the 
most  probable  energy  of  the  ejected  substrate  atoms  when 
oxygen  is  introduced  in  the  A-top  position.   At  present 
there  is  no  apparent  reason  for  such  an  increase.   If  this 
were  a  model  artifact  one  would  expect  the  effect  to  carry 
over,  to  some  degree,  to  crystals  with  adatoms  in  the  bridge 
positions.   Instead,  the  average  and  most  probable  energies 
for  the  bridge  and  three-fold  sites  are  essentially  the 
same. 
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Fig.  15.  Ti-Ti  solid  phase  potential  function  and  associated  parameters 
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Fig.  16.  V-V  solid  phase  potential  function  and  associated  parameters 
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Fig.  17.  Nb-Nb  solid  phase  potential  function  and  associated  parameters 
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Fig.  18.  Ti-0  solid  phase  potential  functionand  associated  parameters 
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APPENDIX  C 
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INITIAL  LOCATIONS  OF  EJECTED  ATOMS 


2. 


Fig.  74.  Origin  of  sputtered  atoms,  clean  Ti 
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Fig.   75.     Origin  of  sputtered  atoms,  Ti   P(2x2)  A-Top 
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Fig.  78.  Origin  of  sputtered  atoms^  Ti  C(2x2)  bridge 
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Fig.   78   (cont.) 
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Fig.  80.  Origin  of  sputtered  atoms,  Ti  C(2x2)TFA 
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Fig.   80  (cont.) 
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Fig.   81    (cont.) 
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Fig.  83.  Origin  of  sputtered  atoms,  Ti  C(2x2)TFB,  0.5  keV 
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Fig.  85.  Origin  of  sputtered  atoms,  V  and  Nb  clean  surface. 
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Fig.  85  (cont.). 
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Fig.  85  (cont.).  Dashed  lines  indicate  second  layer  atoms. 
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Fig.  86.  Origin  of  sputtered  atoms,  V  and  Nb  P(2x2)  A-top. 
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Fig.   86  (cont.). 
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Fig.   87.     Origin  of  sputtered  atoms,  V  and  Nb  C(2x2)  A-top 
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Fig.   87   (cont.). 
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APPENDIX  G 
SPOT  PATTERNS 
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APPENDIX  H 
EJECTION  TIME  VERSUS  ENERGY  CHARTS 
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APPENDIX  I 
SURFACE  LAYER  EJECTION  FREQUENCY  PLOTS 
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Fig.  121.  Relative  surface  layer  ejection  frequency, 
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Fig.  122.  Relative  surface  layer  ejection  frequency. 
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Fig.  123.  Relative  surface  layer  ejection  frequency 
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Fig.  125.  Relative  surface  layer  ejection  frequency 
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EXPERIMENTAL    RESULTS 


Work-function  dependence  of  the  vicld  of  Ti4. 
TiO*.  TiOj*.  and  TiO,*  at  high  oxygen  coverage     The 
Auger  ratio  (or  |0|/|Ti|  is  1.24  tsee  text)     1  he  primary  ion 
was  500-eV  Ne* 


Work-function  dependence  of  the  vield  of  O". 
Tr(>~.  TiOj".  and  TK.),"  at  high  oxygen  coverage     The 
Auger  ratio  lor  lOl/lTi)  is  1.24  tsee  text)     The  exponential 
behavior  was  similar  at  all  other  oxygen  coverages    The  pri 
mary  ion  was  500-eV  Ne* 
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(b»  TiO' 

(c)  Ti02" 

(d)  Ti03" 
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Fig.   127.      (From  ref.    [22]    ). 
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Work  function  dependence  of  the  O", 
NbCr.  NbO,-.  «nd  NbOr  yields.  The  O,  exposure 
was  4  L.   The  primary  ion  was  500  -eV  Ar* . 
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